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^^ ■ The incorporation of priors [ ] in the Optimal Uncertainty Quantification (OUQ) 

framework [35] reveals extreme brittleness in Bayesian inference; very mild assumptions 
O ' are sufficient to demonstrate that, given a set of priors, conditioning on observations can 

■^ ■ produce results as bad as possible, regardless of the sample size. Although it is known 

i*^ , that Bayesian estimators may be inconsistent, if the probability space is non-finite or if 

psj I the model is misspecified, this brittleness result shows that, contrary to popular belief, 

there is no such thing as a "close enough" model in Bayesian inference; a model may 
^ . share an arbitrarily large number of finite-dimensional marginals with, or be arbitrar- 

C/^ i ily close (in Prokhorov or total variation metrics) to, the data-generating distribution 

r~| ' and still make the largest possible prediction error after conditioning on an arbitrarily 

jrt ■ large number of samples. The initial purpose of this paper is to unwrap this brittleness 

mechanism by providing (i) a quantitative version of the Brittleness Theorem of [ i-l] 
(ii) a detailed and comprehensive analysis of its application to the revealing example of 
estimating the mean of a random variable on the unit interval [0, 1] using priors that 
^ , exactly capture the distribution of an arbitrarily large numbers of Hausdorff moments. 

^2 ! However, in doing so, we also found that the free parameter associated with Markov and 

^— s ■ Krem's canonical representations of truncated Hausdorfi^ moments generates reproducing 

t~^ . kernel identities corresponding to reproducing kernel Hilbert spaces of polynomials. Fur- 

t:;^ I thermore, these reproducing identities lead to biorthogonal systems of Selberg integral 

^^ ■ formulas. Therefore, whereas a Selberg integral formula was initially used to compute 

the volume of the Hausdorff moment space, it appears that a refined analysis of the in- 
tegral geometry of the Hausdorff moment space, in the service of proving the Brittleness 
of Bayesian Inference, also reveals a new family of Selberg integral formulas. 
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1 Introduction 

Optimal Uncertainty Quantification (OUQ) [35] provides a framework for the compu- 
tation of optimal bounds on quantities of interest -given a set of available information 
and specified assumptions. Although this framework is neither frequentist nor Bayesian, 
in that it is simply expressed in terms of optimization over measures and functions, a 
natural question arises; what happens when we introduce priors into OUQ? In Owhadi 
et al. [34], the incorporation of priors into OUQ is developed where, in addition to the 
specification of an information/assumption set, a probability measure or a family of prob- 
ability measures is specified on that set. In addition to the development of the necessary 
measure theoretical foundations in the context of Polish spaces, so that simultaneously 
prior measures can be put on subsets of a product space of functions and measures and 
important quantities of interest are measurable, the reduction theory for optimization 
problems over measures on product spaces of measures and functions is also established, 
thus allowing for the computation of optimal bounds on prior and posterior values and 
the analysis of the consequences of conditioning on observed data. 

However, the completion of this program reveals Brittleness theorems [ = , Thm. 5.12, 
Thm. 6.4, Thm. 6.10] for Bayesian Inference -very mild assumptions are sufficient to 
demonstrate that, given a set of priors, conditioning on observations can produce results 
as bad as possible, regardless of the sample size. 

Although it is known from the results of Diaconis and Freedman that the Bayesian 
method may fail to converge or may converge towards the wrong solution (i.e., be incon- 
sistent) if the underlying probability mechanism allows an infinite number of possible 
outcomes [17] and that in these non-finite-probability-space situations, this lack of con- 
vergence (commonly referred to as Bayesian inconsistency) is the rule rather than the 
exception [18], it is also known, from the Bernstein- Von Mises Theorem [10, 47] (see 
also LeCam [ ]), that consistency (convergence upon observation of sample data) does 
indeed hold, under some regularity conditions, if the data-generating distribution of the 
sample data belongs to the finite dimensional family of distributions parameterized by 
the model. Furthermore, although it is also known that this convergence may fail under 
model misspecification [49, 23, 36, 1, 2, 28, 31, 24] (i.e. when the data-generating distri- 
bution does not belong to the family of distributions parameterized by the model), it is 
also a popular belief that a "close enough" model should have good convergence proper- 
ties: see e.g. [19, 40, 20]. The Brittleness theorems [: :, Thm. 5.12, Thm. 6.4, Thm. 6.10] 
show that, on the contrary, there is no such thing as a "close enough" model in Bayesian 
inference in the following sense: Optimal bounds on posterior values obtained from 
models that (i) exactly capture an arbitrarily large finite number of finite-dimensional 
marginals or (ii) or are arbitrarily close (in Prokhorov or total variation metrics) to the 
data-generating distribution show that such models may still make the largest possible 
prediction error after conditioning on an arbitrarily large number of sample data. 

The initial motivation for this paper is to unwrap the mechanism causing brittleness 
in Bayesian inference by providing (i) a quantitative version of the Brittleness Theo- 
rem [34, Thm. 5.12] and (ii) a detailed and comprehensive analysis of its application 



to the informative example of estimating the mean of a random variable on the unit 
interval using priors that exactly capture the distribution of an arbitrary large numbers 
of Hausdorff moments. A thorough analysis of this mechanism is important because 
it implies that, under model misspecification, and without stronger assumptions than 
closeness in the Prokhorov and/or total variation metrics, or in the number of exact 
(generalized) moments, Bayesian inference offers no better guarantee of accuracy than 
arbitrarily picking a value between the essential infimum and supremum of the quantity 
of interest. 

In the example scrutinized in this paper (which corresponds to [34, Ex. 5.15]), the 
probability distribution ^u^ of X is an unknown element of the set of all possible prob- 
ability distributions on [0,1], i.e. /i^ G ^ := 7V4([0, 1]). The set of prior probabil- 
ity distributions vr on /i € ^ (i.e. vr G M.{A)) is defined as the set of priors vr un- 
der which the vector of truncated Hausdorff moments (E^[X], . . . , E^[X"]) is uniformly 
distributed on the truncated Hausdorff moment set M" C W^ defined as the set of 
q = {qi, . . . ,qn) G M" such that there exists a probability measure ^ on [0,1] with 
E^[X*] = Qi with i G {1, . . . ,n}. The computation of optimal bounds on posterior val- 
ues leads us to calculate the Lebesgue volume of subsets of the set M" of truncated 
Hausdorff moments. 

The free parameter associated with Markov and Krem's canonical representations of 
truncated Hausdorff moments (see Section 5) which, along with their principal repre- 
sentations, so handily provides us with the means to prove the quantitative Brittleness 
Theorem 3.3, is found to generate reproducing kernel identities corresponding to re- 
producing kernel Hilbert spaces of polynomials (see sections 6 and 7). Furthermore, 
these reproducing identities lead to biorthogonal systems of Selberg integral formulas 
described in Theorems 8.1, 8.2 and 8.4 (see also Corollary 8.3). 

Therefore, whereas a Selberg integral formula was initially used to compute the 
volume of the Hausdorff moment space, it appears that a refined analysis of the integral 
geometry of the Hausdorff moment space, in the service of proving a Brittleness theorem, 
reveals a new family of Selberg integral formulas (Theorems 8.1, 8.2, 8.4 and Corollary 
8.3). 

On reproducing kernel Hilbert spaces. The first major publication on RKHSs 
is the classic work of Aronszajn [ ], and probably the most famous use of them is 
Bargmann's [7] description of the Fock space of quantum mechanics as a RKHS, the ker- 
nel of which is now known as the Bargmann kernel. Moreover, RKHSs were instrumental 
in the solution of hard and very old problems in Learning Theory, see e.g. Steinwart 
and Christmann [4G] for an excellent treatment. RKHSs are foundational to Splines, see 
e.g. Wahba [ ''^'] and have many applications in Probability and Statistics, see e.g. Berlinet 
and Thomas- Agnan [!•]. Moreover, they are also important in Fluid mechanics, see Be- 
lytschko et al. [ ] and Lui et al. [ ] and Continuum Mechanics, see Lui et al. [32]. 

On Selberg integral formulas. Selberg's integral formulas, n-dimensional general- 
izations of the Euler beta integral, appeared in Norwegian in [ ]. Its importance was 



slow to be appreciated but the recent review "On the importance of the Selberg integral" 
by Forrester and Warnaar [ ] discusses its contributions. From the abstract: 

" We trace its sudden rise to prominence, initiated by a question to Selberg 
from Enrico Bombieri, more than thirty years after its initial publication. 
In quick succession the Selberg integral was used to prove an outstanding 
conjecture in random matrix theory and cases of the Macdonald conjectures. 
It further initiated the study of q-analogues, which in turn enriched the Mac- 
donald conjectures. We review these developments and proceed to exhibit 
the sustained prominence of the Selberg integral as evidenced by its cen- 
tral role in random matrix theory, Calogero-Sutherland quantum many-body 
systems, Knizhnik-Zamolodchikov equations, and multivariable orthogonal 
polynomial theory." 

Regarding the original Selberg integral formulas, according to Forrester and Waardan 
[ , Pg. 3] "For over thirty years the Selberg integral went essentially unnoticed. It was 
used only once- in the special case a = /3 = 1,7 = 2-in a study by S. Karlin and L.S. 
Shapley relating to the volume of a certain moment space, published in 1953." 

2 OUQ with Priors 

To understand OUQ one simply starts with Cebysev [ , Pg. 4] " Given: length, weight, 
position of the centroid and moment of inertia of a material rod with a density varying 
from point to point. It is required to find the most accurate limits for the weight of a 
certain segment of this rod." According to Krem [29], although Cebysev did solve this 
problem, it was his student Markov who supplied the proof in his thesis. See Krem 
[29] for an account of the history of this subject along with substantial contributions by 
Krem. We take this mindset and apply it to more complex problems, extending the base 
space to functions and measures, and, instead of developing sophisticated mathematical 
solutions, develop optimization problems and reductions, so that their solution may be 
implemented on a computer, as in Bertsimas and Popescu's [ ] convex optimization 
approach to Cebysev inequalities, and the Decision Analysis framework of Smith [' ']. 
In addition to the determination of optimal bounds as a function of available information 
and assumptions, the OUQ methodology has the substantial benefit of demanding that 
different components of an organization work together to come up with information and 
assumptions that, together, they believe in. 

Let Af be a topological space and let Q C J-'{X) be a subset of the real-valued 
measurable functions J^{X) on X. Let A be an arbitrary subset of t/ x Ai(X), and 
let ^: Q X M(X) — )■ M be a function producing a quantity of interest. In the context 
of uncertainty quantification one is interested in estimating $(/',^'), where (/',/i') G 
G X A4{X) corresponds to an unknown reality. If A represents all that is known about 
(/^, //^) (in the sense that (/^, /u^) S A and that any (/, /j,) & A could, a priori, be (/^, /Lt^) 



given the available information) then [>-■)] shows that the quantities 

U{A) := sup $(/,m) (2.1) 

C{A) := inf $(/,m) (2.2) 

determine the inequality 

C{A)<Hf\fi^)<U{A) (2.3) 

to be optimal with respect to the available information (i.e. (/',/i') G .4) as follows: 
First, it is simple to see that the inequality (2.3) follows from {f^,fi^) € A. Moreover, 
for any e > there exists a (/, fi) £ A such that 

U{A)-e<m,fi)<U{A). 

Consequently since all that we know about (/^,/i^) is that (/^,//^) G A, it follows that 
the upper bound ^{f',fi^) < U{A) is the best obtainable given that information. The 
lower bound is clearly optimal in the same sense. 

A classical example of a quantity of interest is the validation and certification quantity 
$(/, ;u) := fj,[f > a] where a is a safety margin. In the certification context one is 
interested in showing that /i^[/^ > a] < e where e is a safety certification threshold (i.e. 
the maximum acceptable /iT-probability of the system f' exceeding the safety margin 
a). lil{{A) < e, then the system associated with {f^,fi^) is safe (given the information 
represented by A). If C{A) > e, then the system associated with {p ^fi^) is unsafe. If 
^{A) < e < U{A), then the safety of the system cannot be decided without making 
further assumptions or gathering further information. 

Although the OUQ optimization problems (2.1) and (2.2) are extremely large, we 
have shown in [■'; ] that an important subclass enjoys significant and practical finite- 
dimensional reduction properties. In particular, for assumption sets corresponding to 
linear inequality constraints on generalized moments, the search can be reduced to one 
over probability measures that are products of finite convex combinations of Dirac masses 
with explicit upper bounds on the number of Dirac masses. 

To incorporate priors, we define a prior vr to be a probability measure vr G ^A{A), 
and define the value $(7r) of vr through the extended quantity of interest ^ : A4{A) —?■ M 
defined by 

$(7r) :=E^[$], TTeM{A). 

We will ignore all issues of measurability and refer those interested in this nontrivial and 
not uninteresting topic to [■") 1], which also includes a discussion of more general extended 
quantities of interest. We call the value ]E7r[^] the prior value, and for a family of priors 
n C A4{A) we note that the values 

U{U) := supE^[$] (2.4) 

Tren 

£(n) := infE^r^l (2.5) 



form a natural generalization of the notations U{A) and C{A). Moreover, in the same 
way that U{A) and C{A) are optimal upper and lower bounds on ^{f^fi^) given the 
information that {f^,fi^) € A, W(n) and C{Il) are optimal upper and lower bounds on 
Ett [$] given the information that vr G H. 

For conditioning on sample data in an observation space T>, we begin by defining a 
data map 

B:A^M{V) 

which specifies that B(/, ^) G A4{V) generates the data when the truth is (/, /x) € A. 
Then, given a prior vr G A4{A), we define a probability measure 

7rQB£M{AxV) 

through 

7r0B[^ X B] =E(;.^)^4lA(/,;u)B(/,//)[i?]], A G e(^), B G i3(P) , 

where lyi is the indicator function of the set A: 

We denote the resulting D-marginal by vr • D G 7W (D) which satisfies 

Given an observation d G V, to simultaneously avoid the ill-definedness of regular 
conditional probabilities and incorporate uncertainty in the observation process, we con- 
sider conditioning on an open subset B containing d such that vr • D[B] > 0, where for 
a discussion of the naturality of this positivity condition we refer to [•i-i]- To simplify 
notation, we henceforth drop the notational dependence of the set B on the point d. 
The conditional expectation, given a prior vr and data map D, conditioned on a subset 
B G B{V) such that vr • B[B] > 0, is 

To represent uncertainty regarding the data generating process, instead of a single 
data map ID) : A ^ Ai{'D), we instead specify a set 

D = {B:A^M{V)} 

of data maps and represent our assumptions regarding the data with the statement 
D G 2). Therefore, having specified a set 11 of priors, and a set 1) of data maps, for an 
open subset B CV, we define the set of all possible resulting product measures to be 

n 05 D := |vr D : vr G n, B G D, (vr • D)[S] > o| . 



The notations U(Il) and C{Il) of (2.4) and (2.5) extend naturally to these conditional 
expectations as 

U{UQb^):= sup E^Qo[^\B] 

where we note that, just as for U{A) and C{A) and U{IV) and C{IV), U{Il Qb 2)) and 
/3(n0BS) are optimal upper and lower bounds on the posterior value 'E-kqb [^\B] , given 
the assumptions that vr G 11, B G D, and vr • B(i?) > 0. 

We are now prepared to discuss the brittleness theorems of the next section. Indeed, 
it is easy to see that 

C{A) < C{U) < U{U) < U{A) (2.6) 

and 

C{A) < C{Il Qb D) < U{n Qb 3) < U{A) . 

What the Brittleness Theorem 3.1 will show is that, under mild conditions, regardless 
of where the values C{Il) and U(Il) lie in (2.6) we have 

C{U Qb 3) « C{A) and U{U Qb D) « U{A) , 

that is, conditioning on the observed data, one can obtain any value between C{A) and 
1{{A) for the posterior value E7r0D[*l'|^] for some admissible prior vr € 11 and data map 

BgD. 

3 Quantification of Bayesian Brittleness 

To express the following brittleness theorems, we require, as in [-j-i], the data space T? 
to be metrizable. We select a consistent metric, and for a data point d £ T>, let Bs{d) 
denote the open ball of metric radius 5 about d. To keep the notation simple we omit 
reference to the base point d and denote this family of open balls about d by {Bs, 6 > 0}, 
where Bq = 0. The following generalization of the Brittleness Theorem [34, Thm. 5.12] 
allows a weakening of its assumptions while approximately obtaining its conclusion. 

Theorem 3.1. For a metrizable topological space X , consider a topologized subset Q C 
J-{X) and the space of probability measures M.{X) equipped with the weak star topology. 
Let A C Q X Ai(X) be Suslin, Q separable metrizable, and ^ : A ^ Q Borel measurable. 
Moreover, let £1 C M{Q) be such that supp Q C ^(^),Q G O., and let r > 0. Suppose 
there exists some Q G O, B G S and a continuous monotonically increasing function 
/i : M+ ^ M with /i(0) = such that 

q({q: inf 0{f,fi)[Bs]<T})>l-h{S), 6>0. (3.1) 

Fix 6 > 0. If e >0, e' > and 5' > are three real numbers such that 



sup ^(/,/u)> sup $(/,/x)-y|) >e' 



(3.2) 



and 

then we have 

U{A) - 25' < U{^-\£1) Qbs S) < l^i^) • (3.4) 

If, for r = 0, there exists a 5* > such that for all 6' > there exists some 
Q G Q, B G 3 which satisfies (3.1) with a function h such that h{5) = 0,5 < 5*, and 
which satisfies (3.2) with e = 0, then we recover the conditions and the assertion of the 
Brittleness Theorem [34, Thm. 5.12] for 65,6 < 6*. 

Remark 3.2. The proof of Theorem 3.1 also leads to the following result. For a metriz- 
able topological space X, consider a topologized subset Q C J^{'^) and the space of 
probability measures M.{X) equipped with the weak star topology. Let A C G x A4{X) 
be Suslin, Q separable metrizable, and ^ : A ^ Q Borel measurable. Moreover, let 
Q C M{Q) he such that supp Q C ^(^), Q € 12. It holds true that for 5 > 

U{A) - v{6) < U{^-\£1) Qb, 2)) < U{A) . (3.5) 

where the function v is defined by 



v{S) :=2inf <^(5' >0 



S' > (UiA) - CiA)) inf 

^ ^ (Qe£!,]D)GD,£>0,T>0 



(3.6) 



1 - Q({g : inf(^,^)g^-i(,) ]D)(/,^)[i?5] < r} j + T 
eQ(|g : sup(j^^)g^-i(^)^e(j^^)[B,]>e $(/,^) > sup(j^^)g^ $(/,/i) - (5'| 
for (5 > 0. 

3.1 Application to a revealing example 

To demonstrate that the assumptions of Theorem 3.1 are mild, we now use it to extend 
the Brittleness result of [34, Ex. 5.15] to a simple but informative example. Here one 
is interested in estimating the mean of a random variable X with unknown distribution 
on the unit interval / := [0, 1]. Since our quantity of interest is E t[^], where p) is an 
unknown distribution on /, in the notations of Section 2, we have X := I (since X is a 
random variable on /), Q consists only of the identity function (this example does not 
involve unknown functions of X), A := J^{I) (the set of possible/admissible candidates 
for /x''" is the set of all probability distributions fi on /), $(/x) := Et„.^[t] (our quantity of 
interest is the mean of the random variable X), Q := M" and the map ^ : A4(/) — > M" is 
the map to the truncated Hausdorff moments ^(/i) := (Ej^^[i*]) ._ (our set of prior 

distributions is defined by constraining the distribution of the first n Hausdorff moments 
in M"-, for some fixed n). Furthermore Q is the uniform Borel measure on M"- restricted to 
the Hausdorff moment space M" := ^(A4{I)) and then normalized to be a probability 
measure, that is H C M[A4{I)) is the set of prior distributions on ^ = A4{I) such 



that ^/x € A4{M"^) is uniformly distributed on the the space M"" of first n Hausdorff 
moments. 

The Brittleness Theorem 3.1 imphes (see [•) I, Ex. 5.15]) that if we observe (condition 
on) k independent samples from X, i.e. D := / and B // := fj. ® ■ ■ ■ ® fi (A;-fold 
tensorization) and B is the A:-fold product of small enough balls centered on the data 
then CiJlQB^) ^ C{A) and U{UQb Tl) ^ U{A). In other words, although the set of 
prior values of IE/i[-^] is the single point {^}, the optimal bounds on the posterior values 
of E^ [X] are zero and one irrespective of the number n of constraints on marginals and 
the number k of observed samples if the data is observed with sufficient precision. 

The following theorem provides a rigorous and quantitative statement and proof of 
this implication for k = 1. Although, for the sake of conciseness and clarity our analysis 
is provided in the k = 1 case, it generalizes to the situation where k is arbitrary. Indeed, 
although counterintuitive, one can show that brittleness for the single sample case is 
more difficult to obtain than for multiple samples. Since our main objective here is to 
unwrap and scrutinize the mechanism causing the brittleness of Bayesian inference, we 
therefore chose to keep the presentation and our example as clear, concise, and simple 
as possible to illustrate the generic and pervasive nature of this brittleness. 

Therefore, we will now (i) consider the case of a single data point, i.e., k = 1, V := I, 
and B^/i := fi (ii) use Theorem 3.1 to provide quantitative bounds on ^(11 0^ T>) as a 
function n of the number of marginal constraints defining the set of priors (iii) scrutinize 
the brittleness causing mechanism through the proof of the following theorem. 

Theorem 3.3. Let A := M{I), <^{fi) = Et^f,[t], V = I, and ^ : M{I) -^ W denote 
the map to the truncated Hausdorff moments ^{n) = {^tr~^iAt^])^^i „• Furthermore, let 
Q denote the uniform Borel measure on M"" restricted to the Hausdorff moment space 
M" := ^(7W(/)) and then normalized to be a probability measure. Suppose that Q € Q 
and B^ G 2) . Then for 5 > we have 

2n5 1 
1 - 4e( ) 2"+i < U{^-^{0.) Qbs S) < 1 (3.7) 

Remark 3.4. Alternatively, Theorem 3.3 asserts that for positive 6, 6' satisfying 

6 < l-{6'r^\2e)-'- 

we have 

1-26' <U{^-\£1)Qbs^) <l. (3-8) 

4 Volume Inequalities on the Hausdorff Moment Space 

Karlin and Shapley ['^r;, Thm. 15.2] (see also [27, Thm. 6.2]) computed the volume of the 
space of truncated Hausdorff moments M" of probability measures on the unit interval 
to be 

^<^'') = f\^Tm^^ (4-1) 

10 



where T is the Gamma function. To accomphsh this, they used a Markov representation 
of truncated moment points, as described in Krem [2!)] (see also [27, Ch. II]), combined 
with the change of variables formula, followed by the evaluation of a Selberg integral. 

Here we will refine their analysis to obtain volume inequalities on the Hausdorff 
moment space which are used in the application of the Brittleness Theorem 3.1 to the 
proof of the Brittleness Theorem 3.3. Of the two main results, it is interesting to note 
that the Mass Supremum Equality uses the canonical representation of moment points 
combined with Markov's Maximal Mass Theorem [29, Thm. 2.1] (see also [27, Thm. 4.1]) 
to change the " Inequality" to " Equality" , whereas the Mass Infimum Inequality instead 
uses the principal representation, as in Karlin and Shapley's proof of the volume formula 
(4.1). All this terminology will be defined in the following Section 5 and comes from 
Karlin and Studden [27]. This section will simply state the volume inequalities that we 
need for Theorem 3.3. 

To proceed, let us now fix terminology. Let / := [0, 1], and let 'P{I) be the set of 
Borel measures on / and M{I) C V{I) be the set of probability measures. Throughout 
we will assume the weak star topology for these measures. For the system of functions 

Ui{t) := f,tel,i = 0, ..,n 

the Hausdorff moments of a measure /i € 7^(/) is defined as the vector q € M""^^ with 
coordinates Qi = E^[uj] = Ej^^[t*]. It is well known (see e.g. [1, Cor. 15.7]) that the map 



* : Vil) 



nn+l 



defined by ^(/x) := (Et^^[t'],i = 0,..,nj is affine and continuous. Furthermore, let the 
Hausdorff moment space Al""^^ C M"^^ be the image Al"'*'^ := ^7^(7) of the measures, 
and let M" defined by A4'^~^^ = (1, M") be moments of a probability measures omitting 
the zero-th moment. Equivalently, let Pi : M x M" — )• M" denote the projection mapping 
(xo,xi,... ,x„) onto {xi,...,Xn) and let ^i =: Pi^. Then M" = ^iM{I). We will 
abuse notation by letting ^ also denote the mapping ^i restricted to the first-to-n-th 
order moments of the probability measures 

^ :M{I) ^M" CM". 

It follows from continuity that the moment set 7V4"+^ is a closed convex cone and 
M" is a compact convex set. Moreover, one can show that 

/nt(X"+i) n (1, M") = (1, Int{W)) , (4.2) 

see e.g. [37, Cor. 6.5.1], so that a point q is interior to M" if and only if (1,(/) is interior 
to A^"'"'"^. Let Vol be the usual n-dimensional volume measure. Then, since M" is 
convex, by [12, Lem. 1.8.1] 

Vol{lnt{M'^)) = Vol{M'') . (4.3) 

Our first result is the Mass Supremum Equality. 

11 



Lemma 4.1. Let t^ ^ I, < e < 1, and consider the set M" C M" defined by 

M^ := {g e M" : 3^ € ^'^q : fi{{t,}) > £.} 

Then we have 

Vol{M^) = (1 - e)'"yo/(M") 

Remark 4.2. Lemma 4.1 is valid for any system u, : / ^ M, z = l,..,n of moment 
functions which form a T-system per [27]. 

The second is the Mass Infimum Inequality. 

Lemma 4.3. Let t^ ^ I, 5 > Q, and consider the set M^ C M" defined by 

M^ := {g € M" : 3^ e ^'^q : f^{Bs{t,)) = 0.} 

Then we have 



Vol{M^) >[1- 5(2e)2"j yo/(M") . 
The third is the Mass of First Moment Inequality. 
Lemma 4.4. Let < 6 < g- Then we have 

^ ^^^^ - Vol{M^) -^ 

5 Integral Geometry of the Markov-Krem Representations 

Here we will describe the Markov-Krein representations of truncated moments and begin 
the development of their integral geometry. The history of this subject begins with 
Cebysev and his student Markov's thesis, followed by work by Krem and others, where 
in [29] one can find, not only an historical sketch, but substantial contributions by Krem. 
Indeed, it is clear from Karlin and Studden [27] that this subject owes a lot to Krem. 
Consequently, we refer to the (principal and canonical) representations that we use as 
Markov-Krem representations. It can be argued that the appropriate name should be 
Cebysev-Markov-Krein representations but this name is too long and so we implicitly 
give credit to Cebysev. 

Now, following Karlin and Studden [?", Chapters II & IV], we describe the Markov- 
Krem representations and determine their Jacobian determinants. We finish this section 
by setting up the change of variables approach, in preparation for both the proofs of the 
volume inequalities of Section 4 and all that follows. To wit, we define the index i{t) of 
a strictly increasing set t of points < ti < t2 <■■■< tjsr < 1 hy counting the interior 
points with weight 1 and boundary points with weight 2- For a point q G M" we say 
that a measure /U G A^(/) is a representing measure for q if ^(/x) = q and it is a weighted 
sum of Dirac masses 



N 

z2 ^j"*j ' ^j' 



H = \^X-5^ Xj>0,j = 1,..,N 
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for a strictly increasing set of points < ti < ^2 < ■ ■ ■ < ^A^ ^ 1- In that case, we have 

the formula 

N 

i=i 
The index i(^) of such a representing measure is defined to be the index i{t) of its set 
of support points. 

A representing measure fi is called principal iii{ii) = ^^^^ and canonical if i{^) < ^^^• 
For q E Int[M^), [ , Thm. 2.1] asserts that i(/i) > ^^^ for any representing measure /x 
for q. A principal or canonical representation is called "lower" if its set of support points 
does not include the righthand endpoint 1 and "upper" if it does. The following two 
results will be our main tools. The first is the principal representation, see [27, Cor. 3.1]. 

Theorem 5.1. Every point q € Int{M'^) has a unique upper and lower principal repre- 
sentation. 

The second is the canonical representation which allows the specification of a pre- 
determined point t* G / in the support of the representing measure, see [27, Thm. 3.1] 
combined with [27, Cor. 3.2] and [ , Cor. 3.1]. 

Theorem 5.2. For t^, € (0, 1), every point q G Int{M^) has a unique canonical repre- 
sentation whose support contains t*. When t^, = or 1, there exists a unique principal 
representation whose support contains t*. 

What Theorem 5.2 doesn't make clear is if the canonical representations converge 
to these principal representation as t* tends to or 1. They indeed do as we will see. 
Let us define some notation that we will use henceforward. We consider two coordinate 
representations of the interior of the regular unit simplex. In particular, let 

T^ = {(ti, .., tjv) : < ti < t2 < • • • < tjv < 1} 
denote the set of strictly increasing sequences of length A^ in the interior to / and 

N 

A^ = {(Ai, .., Xn) : A, > 0, j = 1, .., N, J2 ^j < 1} • 

i=i 
denote the interior to the positive orthant restricted to A ■ 1 < 1. Sometimes it will be 
convenient to abuse this notation and shift indices so that 

N-l 

A^ = {(Ao, .., A^_i) : A,- > 0, j = 0, .., A^ - h Yl ^J < ^i ' 

i=o 

We will often use the fact that I^ can be described by A^! copies of T^ corresponding 
to permuting the sequence. 

We use the notation t for a vector with coordinates tj and similarly A for a vector 
with coordinates Aj. We use the superscripts p for "principal" and c for "canonical", 
the subscripts o for "odd", e for "even", / for "lower", and u for "upper". Finally, we 
purposefully ignore multiples of ±1 in all our determinant calculations. With proper 
caution, this causes no harm since at the end of the day we take the absolute value. 
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5.1 Principal Representations 

Theorem 5.1 asserts that each q G Int{M''^) has a unique upper and lower principal 
representation. We now define these representations as maps and compute their Jacobian 
determinants. We state these propositions without proof, since these proofs are very 
similar to those for the canonical representations of Propositions and 5.5 and 5.6. 

First consider the odd case when n = 2m — 1. Then since ^^^ = m is an integer, 
it follows that the support of any lower principal representation contains neither end- 
point and the support of any upper principal representation contains both endpoints. 
Consequently, Theorem 5.1 implies that each point in Int(M^™'^^) has a unique lower 
principal representation of the form 

m m 

^ = ^ Xj6t^ , Xj>0, j = l,..,m, ^ Xj = 1 (5.1) 

i=i i=i 

where < ti < ^2 < ■ • • < i»n < 1- Therefore, consider the bijection 



bli : A™-^ X T*" ^ Int{M^"^~^) 



defined by 



m—l in—l 

\ \ s 

in 



rra— 1 m—l 

^^A,t;. + (l-EA,)C)^J^"^ (5.2) 



^ 2m~l 
2m-l\ 



It also follows that each point in Int{M ) has a unique upper principal representation 
of the form 

m—l m 

fj, = Xo5q + E Xj5tj + XmSi, Xj > 0, j = 0, .., m, ^ Aj = 1 (5.3) 

j=l j=0 

where < ti < ^2 < ■ ■ ■ < Un-i < 1- Therefore, consider the bijection 

(/)P„ : A*" X r'"-^ ^ /nt(M2'"-i) 
defined by 

m—l m—l 

ctP,^{X,t) = m[x,6o+Y^Xj6t^+{l-Y,^,)h) 

3=1 i=0 

m—l m—l „ _, 

3=1 j=0 
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For an increasing sequence tj < ij+i let 

A{t):=llitk-tj) (5.5) 

j<k 

denote the Vandermonde determinant (see e.g. [20, Pg. 400]) of the matrix with entries 
[t*], j = 1, ..,N, i = 0, ..,N — 1. which we write as A^v to emphasize the dimension of 
t. We will also use the same formula for non-increasing sequences when we eventually 
take the absolute value. 

Proposition 5.3. When n = 2m — 1, the Jacobian determinants are 

ra—l m—1 

\det{dcPl,)\{X,t) = J^,{t){l-J2Xj)l[X, 

3=1 3=1 



where 



m—l 

\det{dcpP^)\{X,t) = JP^it) II X, 
J^,{t) := Aiit) 

m— 1 

Note that although each term appears to have the same multiplier n^=i ^j^ i^ ^^^ 
lower case this multiplier is the full product in on A"^~^ and in the upper case it is only 
a partial product on A™, that is, it is missing the Aq term. Finally, let us observe the 
symmetries under the reflection t >—^ 1 — t: 

Jl(l-t) = JLit)- (5.6) 

Now consider the even case when n = 2m. Since ^^^ = m + g is an integer plus -^ it 

follows that the support of any lower principal representation contains the left endpoint 

but not the right and any upper principal representation contains the right endpoint but 

not the left. Let us first consider the lower representation. Theorem 5.1 implies that 

every point in the interior Int{M'^"^) has a unique lower principal representation of the 

form 

mm m 

^ = ^ Xj5t^ + (1 - X] ^J^'^O' ^J > 0' .? = 1' ••' "*' ^^i <^ 
where < ti < • • • tm < 1. Therefore, consider the bijection 
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defined by 

m m 

<t^l,{\t) = ^(^A,<5i^ + (l-^A,)<5o) (5.7) 

X^^f^Z,- (5-8) 



2m 

2m.\ 



On the other hand, every point in the interior Int{M ) has a unique upper principal 
representation of the form 

mm m 

^^ = Y^ ^j^tj + (1 - Z] ^j)^i' ^i > o> i = i> ••! "^. IZ ^i < ^ 

where < ii < • • • t^ < 1- Therefore, consider the bijection 

(/.P„ : A" X T" ^ Int{M^'"^) 
defined by 

m m 

^euiM) = ^,(j2Xj6t,+{l-J2x,)6i) (5.9) 

m m „ 

2m, 



(Ev; + (i-E^^)),ri- (5-10) 



Proposition 5.4. When n = 2m the Jacobian determinants are 



\det{d4^J\iX,t) = J^,it)l[\, 



i=i 



m 

•J 



\det{d<t>PJ\{X,t) = JP^{t)llX 

i=i 
where 

m 
m 

Here, instead of the reflection t i— )• 1 — t leaving the lower and upper invariant as in 
the odd case (5.6), reflection swaps lower and upper; 

J!i{l - t) = J^lit) . (5.11) 
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5.2 Canonical Representations 

Theorem 5.2 asserts that, when t* G (0, 1), every point in Int{M^) has a unique canonical 
representation whose support contains t*, and when t* G {0, 1}, it has a unique principal 
representation whose support contains t#. Therefore, every point in Int{M'^) has a 
unique representing measure ^ such that 

N N 

/" = E^i'^*.' A, >0,j = l,..,iV, J^A, = 1 (5.12) 

3=1 i=i 

such that the sequence < ii < ^2 < ■ ■ ■ < ^Af < 1 contains t*, where for t* G (0, 1), 
the sequence has index ^^^ or ^^, and when i,, = or 1, the index is ^^^. Now let 
us remove t-t from the list and use the identity 'Yl,j=i Aj = 1 to solve for the weight A^, 
corresponding to t^,. Changing notation from N ^^ N + 1 and relabeling the indices, we 
obtain that 

N N 

'" = Ev*. + (i-E^^)'^*- ^^^'^' (5.13) 

where the resulting sequence 

< ti < t2 < • • • < ^TV < 1 

does not contain t*, and when t* G (0, 1), the removal of this interior point implies that 
the resulting sequence has index ^^^ or ^ and when t^, = or 1, the removal of this 
boundary point implies that the resulting sequence has index ^. 

Consequently, for t^, G (0, 1), to represent Int{M'^) we can split into four domains, 
two corresponding to the two ways of producing index ^^^ and two corresponding the 
two ways of producing index ^. When n is even one of the two index ^^^ configurations 
corresponds to including t = in the sequence and not t = 1 and the other corresponds to 
including t = 1 in the sequence and not t = 0, while one of the two index ^ configurations 
corresponds to not allowing t = or t = 1 and the other corresponds to including both 
t = and t = 1. When n is odd this relationships is reversed. Similarly, when t* G {0, 1}, 
we can can split into two domains corresponding to the two ways of producing index ^. 

However, we can show that the representations of index ^^^ produce zero volume 
and so can be excluded from the integral analysis. To that end, we only need to consider 
the ttf G (0, 1) case. Then let us decompose the set of sequences of index N by their 
endpoint configurations. That is, split such sequences into those which contain but 
not 1, 1 but not 0, and 1, and neither or 1. Some of these components will be empty. 
On any of these endpoint specific subdomains let 

Zc{l,..,iV} 

denote the indices of the interior points, so that in this notation we have 

N 

iex je.X'^ i=i 
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Moreover, for a sequence t, let t denote the sequence of interior points, and define 
Tlf := {t : t € T'^ ,tj 7^ t*, j = 1, .., A^} to be the set of interior points which do not cover 
t^ and consider the map 



defined by 



N 
N 

where we note that the first sum ^..gjAjt*- is over the interior points and the second 
^ gjc Ajt* over the endpoints which are fixed. 

The dimension of the domain A^ x Tlf is clearly N + \Z\. However, one can easily 
show that 

N+ \I\ = 2N, 

so that in the case N = ^^^, it follow that the dimension of this subdomain is A^+ |X[ = 
2N = n — 1 < n. Consequently, the image of this subdomain under the map (j) has zero 
volume in M". Since the domain corresponding to index ^^^ is a disjoint union of two 
such subdomains, the assertion is proved. Moreover, the subset consisting of sequences 
which cover t^, also clearly has zero volume, so the constraint that the sequences not 
cover t* can also be removed. 

In conclusion, we can represent the volume VollM'"") using the the representation 

TV N 

J^A,-4 + (l- J^A,)tt, i = l,..,n, (5.14) 

defined on two subdomains corresponding to the two ways that the sequence 

< ti < i2 < • • • < iTV < 1 

can have index ^. That is, when n is even, one subdomain corresponds to not allowing 
or 1 and the other to including both and 1. When n is odd, one subdomain corresponds 
to including and not 1 and the other to including 1 and not 0. 

We now compute the Jacobian determinants. First consider the odd case, n = 2m— 1. 
Then, sequences of index ^ = m — ^ split into the lower and upper sequences 

(0 = ti<t2< ■■■ <t^<l, 
\o <ti <t2 < ■■■ <tm = l- 

Define the lower representation 

(Pll : A™ X r™-i ^ IntiM^""-^) 
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by 

m— 1 m— 1 

4>'oi{X, t; U) = * (Ao<5o + Y, ^A + (1 - IZ ^j)"^** 

m— 1 m— 1 „ _. 

j=l j=0 

and the upper representation 
by 

m— 1 m— 1 

'/'oJA,t;t,) = *(Ao<5i + J^ A,(5i^. + (1 - J^ A,)(5i.^ 

m— 1 m— 1 „ -. 



Proposition 5.5. When n = 2m — 1, for t* € (0, 1), the Jacobian determinants are 

m—l 



\det{d(f>li)\{X,t;t,) = Jo%t„t) II \j 



\detid<PlJ\{X,t;t,) = Jolit,,t) Yl X 



m— 1 

\ 



where 



Moreover, 



m— 1 m—l 

joi{t*,t) := t* n (*^- - **)' n *i • ^--i(*) 

m—l m—l 

J,'L(t*,t) := (1 - t*) n {t, - t,f J] (1 - t,)2 . Al„i(t) 

^oJi*,*) = ^o1(l-i*,l-i)- (5.17) 



Now consider the even case, n = 2m. Then sequences of index ^ = m spht into the 
lower and upper sequences 

Jo < ti < t2 < ••• < im < 1 

[0 = ti < t2 < ••• <im+l = 1, 
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Therefore, we define the lower representation 

cj)li : A™ X T™ ^ Int{M'^"') 
by 



^li{Kt;U) = ^(Y,\A,+{l-Y.X,)dt 



2m 






and the upper representation 

(/.^„ : A™+i X r™-i ^ IntiM^"") 
by 



m— 1 



m— 1 m „ 



j=l j=0 

Proposition 5.6. When n = 2m, for t^ € (0, f), the Jacobian determinants are 



\detidcPli)\iX,t;t,) = JeUt*,t)llXj 



m— 1 

|det(#L)|(A,t;t,) = J,Ut*,t) H Xj , 

i=i 
where 

m 

J,%t,,t):=ll{t,-t,f-Ai{t). 

m— 1 m— 1 

:7'i„(i*,t) := i*(l - t*) n (*J- - **)' n *|(1 - *^)' • ^™-iW ' 

Finally, observe that if we extend the canonical representations to be defined for t* = 
0, 1 by continuity, we obtain the following relations between the canonical representations 
evaluated at the endpoints and the principal representations: 

JoiiO,t) = 

Jon(0,t) = Jo^uit) 
JeMt) = Jl{t) 
Je«(0,t) = 0, (5.20) 
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l#o;(0,i)l 


EE 


l#L(0,i)l 


= l^€Ji)l 


l#e/(0,i)l 


= l#^,(t)l 


l#en(0,t)| 


= 0, 


JoliX.t) 


= JL{t) 


Jon(l,i) 


= 


Je1(l,i) 


= Jluit) 


JeuiX.t) 


= 0, 


Woi{\M 


= I^'/'S.WI 


l#oJl,i)l 


= 


l#d(l,*)l 


= \d€u{t)\ 


l#eJl,i)l 


= 0. 



(5.21) 



(5.22) 



(5.23) 



5.3 Change of variables integral representations 



In Karlin and Shapley's [2(), Thm. 15.2] proof of the Hausdorff moment volume formula 
(4.1), they used the lower principal representation 0^^ of (5.2) when n is odd and 0^^ 
of (5.7) when n is even combined with the change of variables formula. To develop this 
method so that it can be used for the canonical representations, which are not bijections, 
it is convenient to proceed in some generality. To begin, consider a representation 

where W C M" is open and (^ is a continuously differentiable bijection. Then, since 
is injective, by the change of variables formula for injective differentiable mappings 
whose Jacobian determinant may vanish (see e.g. [ J-">, Thm. 3.13] combined with Sard's 
Theorem [15, Thm. 3.14]), we conclude that 

Voi{4>{w))= [ |#|. 
Jw 

Moreover, since (j) is surjective we have 

(j){W) = Int{M") 



and from (4.3) we have 
so that we conclude 



Vol{M'^) =Vol{lnt{M'^)) 



Vol{M'') = / \d^\. (5.24) 

Jw 
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To compute Voli^AI"^), Karlin and Shapley then evaluated the righthand side by de- 
termining the Jacobian determinant and then evaluating the resulting integral using a 
Selberg integral formula. 

However, more can be done along these lines. Indeed, applying the full change of 
variables formula we obtain 

/= / (/o0)l#l 
<I>{W) Jw 

for any function / : (j){W) — > M that is integrable over 4>{W). In particular, since M" is 
compact, it follows using the same reasoning that was applied above to the case / = 1, 
that for any bounded measurable function / : M" — )• R we have 

/= / {fo^)\d<f>\. (5.25) 

M" Jw 

We now apply this to the component functions q i— ?> qi,i = l,..,n on M" where we 
abuse notation and indicate them by the symbol Qi. It may be profitable to also consider 
nonlinear functions such as q ^ qf but we will not do that here. Then, in this notation, 
q- o (p = (j)^ and (5.25) becomes 



qi= / (p'Wl (5.26) 

M" Jw 

That is, we have an integral representation of the mean Hausdorff moments. 

However, to prove Lemma 4.1, instead of a principal representation, we use a family 
of canonical representations from Section 5.2. In this case, utilizing the conclusion at 
(5.14), the major difference with the previous discussion is that, instead of a single 
bijection, there are two continuously differentiable injections 

^k'-Wk^IntiW), A; = 1,2 

that are volume filling in the sense that 

Vol{lnt{W)) = yo/(0i(W^i) U 02(^2)) 

and 

and, instead of Wk,k = 1,2 being open, there exists open sets 14 C Wk,k = 1,2 such 
that 

Vol{Wk) = Vol{Vk) . 

Then the analysis above can easily be repeated to conclude that 

/ /=/ /o0l|#l|+/ /o02|#2| (5.27) 

Jm" Jwi Jw2 

for any bounded measurable function / : M" -^ M. In particular, we conclude 

0l!d</.il+ / 4l#2|, (5.28) 

M" JWi JW2 

our primary integration identity for the mean Hausdorff moments to be used in the next 
section. 
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6 Mean Hausdorff Moments using the Markov-Krem Rep- 
resentations 

We are now prepared to derive integral representations of the mean truncated Hausdorff 
moments with respect to the uniform measure on M" and show that the canonical rep- 
resentations generate reproducing kernel identities corresponding to reproducing kernel 
Hilbert spaces of n-th degree polynomials. These identities are used in Section 8 to de- 
rive biorthogonal systems of Selberg integral formulas. The mean moments with respect 
to many other Selberg-type densities can also be computed but to keep this presentation 
simple we will not do that here. 

We will use Selberg's result (see e.g. [21]) 

i=o 
for the integrals 



T{a + /3 + {n + j- lh)T{l + 7) 



li It 

5„(a,/3,7):= / \[t^-\l - t,f~^\A{t)\^^ dt , (6.2) 

J I'' j=i 

where Re{a) > 0, Re{(3) > 0, Re{-i) > - min (i, Re{a)/{n - l),Re{(3)/{n - 1)) . 

We begin with the volume calculation and then proceed to higher moments using the 
result of the volume calculation. The main idea of our approach is the following. Recall 
from Section 5.1 that the lower and upper principal representations are each bijections 
with Int{M"^) so that the volume Vol{M'^) can be computed using the change of variables 
result (5.24). For example, when n = 2m — 1, the lower principal representation 0^^ 
defined in (5.2) and the upper principal representation (/>^^ defined in (5.4), along with the 
values of their Jacobian determinants from Proposition 5.3 produce two different integral 
representations for Vol{M"'). Specifically, in the notation for Selberg's formulas (6.1) 
for the integrals (6.2), using the identity /^™^_i (l - Ylf=i^ ^j) llT=i^ '^j^'^ = (2m- 1)! ' 
the lower representation yields 

VoliM^""-^) = [ |det(#^^)| 

jn— 1 m~l „ 

(i-E^.onvA)/ J!, 

1 f' J, ' 



{2m - 1)! Jt. 

(2m-l)!m! /^ ^"'^ ^^'^^ 

i2m-l)lmd ^^-^'^'' 
^ 5m(l,l,2). 



(2m — l)!?7i! 
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On the other hand, using the identity J^^ TYjLi ^jdX = (2m-i)\ ' ^^^ upper representa- 
tion (jj^u yields 



rn—l „ 

n ^j'^A / ^ou 



rn—l 

A" 

j — ^ 

1 {■ 

JP 

^ mi. 



(2m — 1)! Jt"^-i 

(2m-l)!(m-l)! Jjr^-i ^°" ^^'^^ 

{2m-l)\{m-l)\ Jjm-i ^J- ^^ ^' ™^ ^ 
r'S'm-i(3, 3,2) . 



(m-l)!(2m-l)! 
Combining the two results (6.3) and (6.4) we conclude the identity 

r'S'm-i(3, 3,2) = — -S'm(l,l,2) 



[m — 1): m\ 

which is confirmed through direct calculation. 

In the even case, where n = 2m, we use the representation 0^^ defined in (5.7) and 
its Jacobian determinant from Proposition 5.4, along with the identity /^^ TlfLi ^jdX = 
/^^N, , to conclude that 

voiiM^n = —\- f J!i = 7:r^r-M^^ i' 2) ■ (6-5) 

{2m)\m\ Jim {2m)\m\ 

Using same identity, the upper representation ^f^ defined in (5.9) yields 

^°'(^^'""> = » I,, "'" = (S^^^-'i' ^- ^> ■ <"■<" 

Equating the two we conclude that 

5^(1,3,2) = S„,(3, 1,2) 

which is well known from the symmetry of the Selberg formula in its first two arguments, 
and corresponds to the change of variables t i— t- 1 — t. Consequently, we see how two 
different integral representations of the volume Vol{M'^) generate identities. 

However, the canonical representations form a one parameter family of representa- 
tions of Int{M^) and the value Vol{M^) expressed in terms of the resulting one param- 
eter family of integrals produces more interesting results. To see this, consider the odd 
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case n = 2m — 1, and the volume filling pair of representations (/)^j and (j)^^ defined in 
(5.15) and (5.16) with Jacobian determinants evaluated in Proposition 5.5. Apply the 
modified change of variables formula (5.28) in terms of these two representations, and 
the identity /^,„ Hjll ^jd^ = (2m-i)\ ' *° obtain 

Vol{M^^-') = [ |#^,| 

m— 1 „ „ m— 1 



(2m -1)! /r™-i ^^"' ^ ^°"'' 



c 



(2m-l)!(m-l)! Jim-i 
Therefore, showing the parameters, we conclude that for t^, E (0, 1) we have 

Since the identity (6.7) holds for all t* G (0, 1) it generates integral identities. For the 
first, since the integrand is continuous in t* we can set t* = to obtain 

^°'(*^'""' = (2„-l)U-l)! L, (■^"'C' *> + ■^•"(°- '»* • 
but from Propositions 5.5 we have 

Jo;(0,t) = 
and from (5.20) 

JLi^.t) = JL{t) 

so that we obtain 

which we already knew from the volume calculation using the principal representation 
(6.4). However, if we compute the first order differential invariant by differentiating (6.7) 
with respect to t^, at i* = we obtain the first integral formula of Theorem 8.1. 

Now consider the even case n = 2m, and the volume filling pair of representations 
(j)'^i and (/)g^ defined in (5.18) and (5.19) with Jacobian determinants evaluated in Propo- 
sition 5.6. Apply the modified change of variables formula (5.28) in terms of these two 
representations, and the identities /^™. Hjli ^jd^ = (2^ ^^^ /a'"+i Hjl'i ^id^ = (2^ 
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to obtain 



Vol{M^^) = f \dct>li\ + f |#^J 

JA'"xT'" JA'^+ixT™-! 

m „ „ m—1 



n ^j-^^) ^ei+{ n ^^■^^) / ^- 



A-n 



(2m)! y^m "' (2m)! yj.™-! 



{2m)\m\ J im ""' {2m)\{m-iy. Jjm-i 
Showing the parameters, we conclude that for t^, € (0, 1) we have 

VoliM^"") = . \ , / J^i{U,t)dt + — — -^ — / JeUt,,t)dt . 

{2m)\m\Jim {2m)\{m - 1)\ J jm-i 

Setting i* = and using 






from Proposition 5.6 and 

from (5.20), we obtain 

VoliM^"') = j;r^ [ J^i(.0,t)dt + -—^ -[ X^{0,t)dt 

(zmjim,! Jjm. (2mj!(m — Ij! Jjm-i 



1 



(2m)!m! 



J!l(t)dt 



(2m) 
1 






(2m)!m,! 



5m(3,l,2) 



which we ah'ead knew from the volume calculation using the principal representation 
(6.5). However, if we compute the first order differential invariant by differentiating 
(6.8) with respect to t* at t* = we obtain the second integral formula of Theorem 8.1. 
We can now proceed to compute the mean of the moments with respect to the 
uniform measure on M" using the volume identities (6.3), (6.4), (6.5), (6.6) from the 
principal representations and (6.7) and (6.8) from the canonical representations. From 
the identities (5.20), (5. 21), (5. 22), (5.23) connecting the canonical representations at the 
endpoints and the principal representations, it is clear that we can generate the integral 
formula for the mean moments corresponding to all the principal representations except 
(f)^i by doing so using the canonical representations and then evaluating the result at the 
endpoints. Therefore, we move directly to the canonical representations. Let 5o denote 
the indicator function defined hy 6o{i) = l,i = and 6o{i) = otherwise. Using the 
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convention that O'^ := 1, the fohowing proposition utihzes the volume equalities (6.7) and 
(6.8) to simultaneously expresses themselves and the moment equalities generated by the 
canonical representations. For a function (/> : / — ?> R we define the diagonal extension 

S(/. : /^ ^ R by 

N 

(S0)(t):=^0(t,), tei^. 
i=i 
For simple powers, we introduce the notation 



Ef:=^f 



N 

i 
J 



for the power sum and note the important example 

N 

that will be used in the Selberg integral formulas of Theorem 8.1. 
Proposition 6.1. Let n = 2m — 1. Then for all t^, ^ I and i = 0, 1, .., 2m — 1 we have 

t\ 



f Qi- ^Vol{M 



2m-lN 



^o(^) f ^cu .^.. , 1 



- i2m)Km - 1)! L. ^-^'-'^'' + (2..)!(1 - 1)! L-. ^-^'-'^''' 
+ (2^1(1-1)! L.^<-^"'(^-^) + -^-(^-^))^^- 



Let n = 2m. Then for all t^, ^ I and i = 0,1, .., 2m we have 
«^ - ^T^rVoliM'n 

^/2m 2m + 1 

~ (2m + l)!(m-l)!A™-i^^"^*'^'' 

+ 7^ h^^^-^ ^t'Jel{t*,t)dt + j:: — ^ ttt/ ^fj^^{t^,t)dt. 

{2m + ly.ml J j„, "'■ (2m + l)!(m - 1)! J/m-i 

The above technique of comparing two representations of the same volume to gen- 
erate identities we now employ by simply subtracting the integral representations of 
Proposition 6.1 evaluated at t* = from that with arbitrary t^ € I. We now show that 
this procedure produces a reproducing kernel identity on the space of polynomials. 
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7 The Canonical Representations and Reproducing Kernel 
Hilbert Spaces of Polynomials 

The integral representations of Proposition 6.1 show clear signs of the existence of re- 
producing kernel identities of the form 



fix) = j K{x, y)f{y)dy, f e H, x e X 



since, in the odd case, the integrand on the righthand side St* is integrated against a 
kernel J^giit*,t) + J'^^{t^,,t) and produces a multiple of tl plus some terms. Reproducing 
kernel Hilbert spaces are Hilbert spaces of functions such that pointwise evaluation is 
continuous on the Hilbert space. They have remarkable properties, in particular, the 
reproducing kernel identities which can be thought of like an abstract Cauchy integral 
formula from complex analysis. 

Let us present Proposition 6.1 in reproducing kernel form. To that end, define 

nt*,t) ■■= Jou{0,t) - JoUU,t) - J^{U,t), (7.1) 

and note that from (5.20) we have J^oui^^'^) — J^ouit), so that 

m*,t) = JL{t) - Joiit*,t) - JLiu.t) . 

Moreover, observe that the symmetries 

of (5.6) and (5.17) combined with J'^i{0,t) = imply that 

'H{0,t) = 

n{l,t) = 0. (7.2) 

and 

n{i-t^,i-t) = n{u,t). (7.3) 

Let H" denote the space of n-th degree polynomials in one variable with real coeffi- 
cients. 

Theorem 7.1. For all (p G n^m-i y;e have 

4>{U)Vol{M^^-^) = -^ / {^cp){t)n{U,t)dt 

(2m — l)][m — 1)! Jjm-i 

+ r^ tw Tv/ '^^ait*,t)dt (7.4) 

[2m — l)!(?7i — 1)! Jjm-i 

+ To ^^ TTt/ ^oi{t.,t)dt 

[2m — l)\[m — Ij! J/m-i 

28 



and for (j) G n we have 



[2m)\[m — Ij! Jjm-i 



/ J,%U,t)dt - f/^^ / Jeuit*,t)dt . 

J/m (2m)!(m- Ij! J/m-i 



(2m)!m! j/m 

To integrate out the diagonal extension S, for any function (t*,t) i-^ J'(t*,t) we let 
J{t^,s) := / j{t^,{s,t2,..,tN))dt2 ■■■dt]\f 



denote the marginalization to the first component of t. Now for any such function J, 
which is invariant under the symmetric group acting on its second variable, we have 



{i:^){t)J{U,t)dt = N / (l){s)J{U,s)ds 
i^ Ji 

so that we obtain the following corollary to Theorem 7.1. Let us define 

GiU, s) := JMO, s) - MU, s) ~ {m - 1)J^^{U, s) 
and note that 

g{o,s) = o 

but 

g{i,s)^o. 

Let Hq C n" denote the n-th degree polynomials (p € H" which vanish on the boundary 
of /, that is, 0(0) = </)(l) = 0. 

Corollary 7.2. For all (f> G n^™^-'^ we have 

4>{U)Vol{M'^^-') = -^ f^^s)n{U,s)ds 

[2m — lj!(?TT, — 2)] J J 

+ To tw ^J ^ou{U,s)ds (7.5) 

[2m — l)!(?7i — 2)\ J I 

+ To ?S^ l^J Jol{t*,s)ds 

[2m — 1)!(?TT, — 2)\ J I 
and for (/) G n^™" we have 

4>{U)Vol{M^n = ^2^y_^^_^yJ ^4>{s)Q{U.s)ds 

{2m)\{m - ly. J J {2m)\{m - 2)\ J J 
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In particular, for the normalizations 



we have 



yo/(M2'»-i) (2m - l)!(m - 2)! 
yo/(M2™) (2m)!(m - 1)! 



(/,(t,) = (f>{s)n{t,,s)ds, (/.en^'"-! 

</>(**) = [ (f>{s)g{t„s)ds, (j)€Ul'^. 



Let us now restrict our attention to the odd case and let LP' {I) denote the usual 
Lebesgue space corresponding to the uniform Borel measure on /. Then, it is well 
known, see e.g. Saitoh [ , Thm. 1, Pg. 21], that the integral operator 



^ f (Pis)nit^,s)ds, (pGL'^il) 



determines a reproducing kernel Hilbert space structure on its range with reproducing 
kernel 

/C(ri,r2) := / 'H{ri,s)'H{r2,s)ds . 

From the definition (7.1) 

nt*,t) := JS.(0,t) - J.Kt,,*) - JS,(i*,t) 

and the definitions of ^^^ and JTJ^ from Proposition 5.5, it follows that 'H{-,t) E n^"*"^, t G 
/™-i, and therefore it follows from (7.2) that n{-,t) G n^"^"!, t G /'"-i. Consequently, 
by marginalization to % and scalar normalization, we have 

7^(-,s)Gn2™-i, sGl. (7.6) 

Therefore the range of this integral operator is contained in Hq™" . However, it follows 
from Corollary 7.2 that the range is identically Hq™" . Therefore we conclude that 
IIq™" is a reproducing kernel Hilbert space with kernel /C. 

Because of Corollary 7.2, one might be tempted to think that this reproducing ker- 
nel Hilbert space structure corresponds to that which Hq™" inherits as the subspace 
IIq'"" C L'^{I), but this is not the case. Indeed, let P2m-i denote the L'^{I) orthogonal 
projection P2m-i : L'^{I) -^ ^o^~ ^iid consider the kernel 

IC2m-i{ri,r2) := / P2m-i'H{ri, s) ■ P2m-i'H{r2, s)ds (7.7) 

where the projections are acting on the kernels in the second component. Then, since 
this projection makes no difference in the reproducing identities in Corollary 7.2, one 
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can show that K,2m-i is the reproducing kernel associated with Hq™^ C L'^{I) and since 
the latter can be computed in terms of the Legendre polynomials of order 2 (see e.g. [5, 
Sec. 12.5]) using the Christoffel-Darboux formula [15] (see e.g. Simon [42] for a more 
current reference), we conclude an identification of /C2m-i with the Christoffel-Darboux 
formula for the kernel of the Legendre polynomials of order 2. That, is 

^^"^-^^"^'"^)-2(2m-l)(2m)(2m + l) ^7^7^ ^^'^^ 

for (ri,r2) G /^, where Pk are the Legendre polynomials shifted to the interval 



1 d'^(r" — r 



k 



and 

Qfc(r):=(r-r2)P^'(r), rG/ (7.10) 

are the associated Legendre polynomials of order 2 (see e.g. [ ,, Sec. 12.5]). 
Moreover, since 

H{U,s) := J„^JO, s) - J^i{U,s) - J,liU, s), 

and from Proposition 5.5 we have 

^o«(0,0) ^ 

JSuit.,0) > 0, i*G(0,l), 

we find that 

n{U,0) = -j:;^{U,0)<0, t*G(0,l). 

Consequently, for t^ G (0, 1), it fohows that 'H{t^,, •) ^ Hq™""^ and therefore 

K- / /C2m-1 • 

Moreover, from the orthogonal decomposition 
JC{ri,r2) = 'H{ri,s)'H{r2,s)ds 

= / P2m-i'H{ri,s) ■ P2m--iH{r2, s)ds + / P2m-ii-i{ri, s) ■ P2m-iH{r2, s)ds 

= }C2m-i{ri,r2) + / P'L.-i'^ ■ P'L.-i'^ir2,s)ds 
we conclude 
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Theorem 7.3. Let K,2m-i denote the reproducing kernel for the polynomials Hq™'" as 
a subset of L^{I). Then /C2m-i can be expressed by both (7.7) and (7.8). Moreover, 
riQ™" is also a reproducing kernel Hilbert space with kernel K,, and 

K. — /C2m-1 

is a reproducing kernel. 

Remark 7.4. Besides the fact that the kernel /C defining the Hilbert space structure for 
the polynomials Hq"^" is not that of the Legendre polynomials, we do not know if this 
kernel is known, nor do we have an explicit formula for it. However, what this section 
shows is that this kernel and its associated Hilbert space Hq*"" are intimately connected 
with the canonical representations of truncated Hausdorff moments, and therefore might 
be called the Markov-Krem kernel. Moreover, if instead of the uniform measure on the 
moments, a Selberg type density is used, more such reproducing kernels may be revealed. 

8 New Selberg Integral Formulas 

The integral representations of the mean Hausdorff moments of Proposition 6.1 provide 
new integral identities of Selberg type. In the following theorem, we provide the first 
in a sequence corresponding to when n is odd and even. We then show how to use the 
reproducing kernel identities of Theorem 7.1 to generate biorthogonal systems of Selberg 
integral formulas. 

Theorem 8.1. It holds true that 

[ St-i .f[t]il - t.f^tm = ^^(5. 1.2) -^^(3, 3, 2) _ ^^^^^ 

and 

/ i:t-^-Wt]-^tn{t)dt = ^S^^i{f>,^,2). (8.2) 

The identities of Theorem 8.1 follow only from the volume equalities, that is, the 
i = case of Theorem 7.1. The following theorem demonstrates how to use all the 
moment equalities of Theorem 7.1 to generate biorthogonal systems of Selberg integral 
formulas. Let us recall definition (7.1) 

nu.t) = JIM - j,M,t) - joiit.,t) . 

Theorem 8.2. Let n = 2m — 1 and consider the scaled kernel 

y .^ _^ 2 

■ yo/(M2™-i)(2m-l)!(m-l)! ' 

Then, 

n{;t)eUl"'-\ tel"'-^ 
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Moreover, consider a basis {pj,j = 1, ••, 2m — 2} for Hq™ and the resulting expansion 
of'H{-,t) in this basis for each t E I"^~^ ; 

2m-2 

'H{t^,t) = Y^ hj{t)pj{t^), (t*,t) G / X r'-^ . 

Then, {Tipj,j = l,..,2m — 2}, {hj,j = l,..,2m — 2} form an L^(/™^-'^) biorthogonal 
system. That is, 

hjEpk = 5jk, j, k = l, ..,2m - 2 



f Jm—l 

As an immediate corollary, we have 
Corollary 8.3. Let n = 2m — 1 and consider the scaled marginal kernel 

fl .= _^ 2 

' yoZ(M2'»-i)(2m-l)!(m-2)! 

(note the different scaling than Theorem 8.2). Then, 

?^(-,s)Gn2"^-\ s&i. 

Moreover, consider a basis {pj,j = 1, ..,2m — 2} for Hq™^ and the resulting expansion 
of'H{-, s) in this basis for each s S I; 

2m-2 

1-i{t*,s)= "Y hj{s)pj{t^:), (t*,s)E/x/. 

i=i 

Then, {pj,j = l,..,2m — 2}, {hj,j = l,..,2m — 2} form an L'^{I) biorthogonal system. 
That is, 

hjPk = Sjk, j,k = l,..,2m-2. 

The choice of basis for Hq™"" determines the corresponding component functions 
hj,j = l,..,2m — 2 and the integrands hjTip^ in Theorem 8.2. Therefore, the task 
remaining is to select a basis for which the component functions hj can be determined and 
such that the resulting integrals are of interest. When the chosen basis is orthonormal 
with respect to some inner product (•,•), then the coefficients hj in the representation 

2m-2 

'H{t^,t) = Y^ hj{t)pj{t^), {t^,t) el X /™"^ . 
i=i 

of Theorem 8.1 are 

hj{t) = {n{;t),pj) . 

As an example, we now compute these component functions, and therefore deter- 
mine explicit forms for these Selberg integrals, when the basis consists of the associated 
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Legendre polynomials of order 2. To that end, recall the definitions (7.9) and (7.10) of 
the Legendre polynomials and the associated Legendre polynomials of order 2 translated 
to the unit interval /. In addition, recall the j-th symmetric function ej defined as 

ej[t) := y ^ ti^- • • tij 

with Co := 1 and the symmetric functions ej{t, z) restricted to the diagonal t = z 

ej{t,t):= Yl ^h(t)('nii)^ j = 0,..,2m-2. (8.3) 

jl+J2=j 

Theorem 8.4. Consider the basis of IiQ™'~ consisting of the associated Legendre polyno- 
mials Qj,j = 2, .., 2m— 1 of order 2 translated to the unit interval I. For k = 2, .., 2m — 1 
define 

^^ r{j + k + 2)r(j - A; + 1) ' - -^ - 

2m- 1 

hk{t) := ^ {-iy^'^ajke2m.-i-j{t-,t). 

j=k 

Then for j = k mod 2, j,k = 2, .., 2m — 1, we have 



j^^_^ /i,(t)SQ,(t) n 4 • ^t-iit)dt = Vol{M'^-'){2m-mm-l)lj^^^^±^ 



m— 1 

„ . „ , ( Ac -I- ■/ 1 1 



9 Proofs 

9.1 Proof of Theorem 3.1 

We seek to apply the nested reduction theorem ['! ', Thm. 5.10]. The assertion is trivially 
true when C{A) = U{A) so we can assume C{A) < hl{A). Let us first establish that 
the assumptions of the theorem are well defined. To that end, note that [•) i:, Lem. 4.11] 
(which follows from Castaing and Valadier [l-''), Lemma in.39 p. 86], which in turn fol- 
lows from Saint-Beuve's [ . ] extension of Aumann's Selection Theorem to Suslin spaces) 
implies that q — )■ inf(j^^)g^-i(q) D(/, ^)[i?] is universally measurable and hence the con- 
ditions of the theorem are well defined if we extend the definitions in the usual way when 
operating on universally measurable sets and functions. Similarly, since for any A the 
function (/, /i) i-^ (<!>(/, /i) — A)B(/, /x)[i?] is measurable, the function ^ : Q — > R defined 
by 

9{q):= sup ($(/,//)- A)D(/,^)[i35] 
(/,/.) Gvf- 1(g) 

is universally measurable. 
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For the proof of the theorem, let Q €1}, B € 3 satisfy the assumptions, and define 
A := U{A) — 26' . Consider the events 

Q5:={q-.^ inf B(/,/i)[i?5]<T} 

Ue,s = \q ■■ sup ^(/,Ai) > sup $(/,^) - 6'} 

where the assumptions (3.1) and (3.2) become 

Q{Q5)>l-H5), 5>0 

Q{Ue,5) >e', d>0 
Let us denote r' := t(U{A) — C{A)). It is easy to see that 

{0 > -t'} D Qs 

and 

{9 > s5'} D U,,s 

and therefore 



and 



> 1 - h{5) 



>eS'}) > q{Ue,5) 
> e' . 



Since <!>(/, ^) > C{A), {f-,pi) G A, it follows that [6*1 < U{A) — C{A), and so we obtain 

f edQ = f edq+ I 9dQ 

J J{e>e&'} J{9<eS'} 

> e5'q{{e > e6'}) + / 9dQ 

J{e<eS'} 

> e5'q{{0 > e5'}) + / Odq 

J{e<o} 

> e5'q{{e > e6'}) + [ 9dq+[ 0dq 

J{e<-T'} J{-r'<e<o} 

> e5'q{{e > e5'}) - {U{A) - C{A))q{{9 < -r'}) - r'Q({-r' <9<0}) 

> eS'e' - {U{A) - C{A)) h{8) - t' . 

Therefore, for any strictly positive solution 5 > to 

eS'e' 



h{6) + T < 



U{A) - C{A) 
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we have 

Eq[9]= f edQ>0, 

where we recall that the function 6 depends on 6, and therefore trivially 



sup Eg^ 

QeQ,ID)G2) 



sup {m,fi)-X)D{f,fi)[Bs] 
(/,A.)e*-i{g) 



>0. 



The assertion then follows from [ , Thm. 5.10]. 



9.2 Proof of Theorem 3.3 



We will apply the Mass Supremum Equality 4.1, the Mass Infimum Inequality 4.3, and 
the Mass of First Moment Inequality 4.4. To that end, define the events 

5,,5 = {g G M" : 3/i G ^"^g : ^,{Bs) > e} 

Is = {q€ M" : 3^ G ^-\ : f^{Bs) = o} 

FM^, = {gGM":(7iG(l-5',l]} 
First observe that some endpoint conditions have zero mass. For example, 

Q{Se,5) := Q{{q : 3/z G ^''q : i,{Bs) > e}) = Q{{q : 3/x G ^"^g : fi{Bs) > e}) 
and 

Q{FMs>) := Q{{q G M" : gi G (1 - 6', 1]}) = Q{{q G M" : gi G [1 - 6', 1]}) . 
Consequently, the Mass Supremum Equality 4.1 asserts that 

^e,5) = Q{{q ■■ 3/U G ^-'q : fi{Bs) > e}) > (1 - e)^ 



where the right-hand side is independent of 5, the Mass Infimum Inequality 4.3 asserts 
that 

q{Is) > 1 - 6{2ef^ , 

and the Mass of First Moment Inequality 4.4 asserts that 



Q{FMs> 
Define the events 



qeM^:q,e[l-5',l]})>{6'r. 



Qs := {q : inf fi[Bs] = O} 



Ue,5 ■■= \q : sup E^[X] > 1 - 5'} 
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Then since 
we have 
and since 






we have 



Ue,5 = \q: sup E^[X]>l-6'} 

^ /.G*-i(g),M[Bi]>£ ^ 

= {^Z : 91 G (1 - 5', l],3f,G ^-'q : fi[Bs] > e} 

= S.sriFMs' 



{Us,s) = q{Se,5nFMs') 

= l-Q{{S,,snFMyr) 

= i-q{sis)-Q{FM^') 

= Q{Ss,5)+Q{FMs')-l 

= (l-e)"-l + (y)" 

= {6T-ne 



Consequently, if we choose e := ^hrj—, then 



2n ' 



2 ' 
so that the assumptions (3.1) and (3.2), expressed as 

q{Q5)>i-H5), 5>o 

Q{Ue,s)>e', 6>0, 
are satisfied with s' := ^^-^, e := 2n ' ^^'^ ^('^) •" (^(26)^" . We can solve 



U{A)- 


-C{A) 


eS'e' 




isr,, 


{ST 


2n ^ 


2 


(j/-)2n+l 



An 

1 /e/\2n+l/r, \— 2n 



by choosing 6 < -^{^') " (2e 
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9.3 Proof of Lemma 4.1 

It follows from the discussion above Equation (5.14) in Section 5.2 that the integration 
formula (5.27), applied to the indicator function 



of the set 

can be expressed 



M^ := {g G ikf^ : 3/i G ^-\ : n{{t,}) > e} 

Vol{M^) = Y^ Vol{M^ n MW')) , (9.1) 

i=l,2 

where (j)k,k = 1,2 correspond to the two ways that the sequence of interior points T 
can have index ^ in the representation 

defined by 

TV 

where we note that the first sum X^iei'^i^i' ^^ '-'^^^ ^^^ interior points tj G T and the 
second X],gxc Ajt*- over the endpoints which are fixed. In particular, when n is even, one 
subdomain corresponds to not allowing or 1 and the other to including both and 1. 
When n is odd, one subdomain corresponds to including and not 1 and the other to 
including 1 and not 0. For each such representation, we have 

N+\I\ = n. (9.3) 

Consider one of these representations, and let A^ C A , defined by 



N 

Af:={AGA^:j:A,<l-e} 



define the subset 

Af X f^ c A^ X f ^ 

corresponding to those moment points represented with the mass allocated (1 — X]'i=i ^^i) 
to the point t* being greater than or equal to e. It is clear that 

(/>(Af X f^) c AC n 0(A^ X f ^) , 

however, the canonical representations have an extremal characterization by Markov, 
see Krem [:''\ Thm. 2.1] (see also [^'^, Thm. 3.1]), which asserts that 

</.(Af X f^) = M," n (/<(A^ X f ^) . (9.4) 
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It follows from the affine dependence on A in (9.2) that the determinant \det{d(p)\ 
splits into a product 

|(iet(#)(A,t)| = (T[Xj)\det{d(l)){l,i)\ 

so that integration also splits into a product 

VolU{A^ xf^)) = f |(iet(#)| 

JAN^fN 

= / \det{d(f>){X,i)\dXdi 



\Ns^rpN ^^ 



(TlXj)\det{d(l)){l,i)\dXdi 

iex 

Wx^dX I \det{d(t>){l;i)\di 

^_ ^JV+|I| /■ (Y\x\dx[ \det{d(^){l;i)\di 
Jan ^-K-i- ^ Jfiv 



^^ ..X 



J ex 

= (l-e)Vo/(</.(A^xf^)) 

where the last equality follows from (9.3). That is, we have 

yo/(0(Af X f^)) = (1 - e)"yo/(</.(A^ X f^)) . (9.5) 

To combine the results for both representations (pi and (j)2, let us denote the correspond- 
ing domains W^ and W'^ of the form A x T and the corresponding subdomains W^ 
and W^ of the form A^ x T^ . Then, the maximal mass identity (9.4) becomes 

Vol{M,"nUW')) =Vol{cl)iiW^)), i = l,2, 

and (9.5) becomes 

Vol{(t>{W',)) = (1 - ey'Vol{<p{W')) i = l,2. 

Consequently, from the integration formula (9.1), we conclude 

Vol{M^) = ^Vol{M^n(l)i{W')) 

i=l,2 
i=l,2 
i=l,2 

= (i-.)"^yo/(0(w^^)) 

1=1,2 

= (l-e)"Fo/(M") 
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Therefore, the assertion is proved. 



9.4 Proof of Lemma 4.3 

First consider the odd case, n = 2m — 1. We utihze the bijective principal representation 
0^; : A'"-^ X T*" ^ IntiNP"^-^) defined in (5.2) and 

m—l m—1 



|ciei(#^,)|(A,t) = J-^(t)(l-^A,)nA, 



J^M := Ai{t) 



where 

from Proposition 5.3 along with the change of variables formula (5.25). 
Fix t^ G (0, 1) and let 

TP := {(ti, .., t„) G T'- : t, ^ Bs{Q,j = 1, ..,m } 

denote those sequences which have no point a distance less than 6 from t^. It follows 
that 

and therefore 

Fo/^M^"-!) > Fo/^^P (A"-i X Tf )) . (9.6) 

We bound the righthand side from below using the change of variables formula (5.25) as 



l-o/te(A™-ixrr)) = / \det{dcPl,)\ 

„ m—l m—l „ 

7 = 1 7 = 1 5 



and then bounding 



I J^idt = [ Ai{t)dt 



AUt)dt 



ml I T-n 

where 

If := {(ti,..,tm) G I'" : tj ^ B5{U),j = l,..,m}. 

To bound this from below we bound the integral over (/™) from above. To that end, 
let 

I^^ := {(ti,..,tm) G I™ : tj G i?5(i*)}, j = l,..,m. 
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so that 

Therefore, using a union bound and the symmetry of A we have 
/ ^^Ai{t)dt = [ At{t)dt 

m „ 

= m A^{t)dt 
'''s.i 

= m Y\ {tk-tj)^dti---dtm 

■S'l l<j<k<m 

< ml W {tk-tj)'^dti---dtm 

I jm -*■-'- 

*.! 2<j<k<m 

= mVol{Bs) I W {tk- tjfdt2 ■■■dtr. 



2<j<k<rn 

= mVol{Bs) I A'^-iit)dt 

J im-l 

< 2m55™_i (1,1,2) 



and so obtain 

VolUiA"^-^ X T"*)) //- ^fnii)(^t 

> l-2n..%#Ml 
5^(1,1,2) 

Using Selberg's formulas (6.1) we compute 

Q n 1 0^ ■prm-2 r(l+2im3+2i) 

'^m-l(l,l,2j _ llj=0 2r(2(m+i)-2) 



5^(1,1,2) T-rm-l r{l+2jfr{3+2j) 

^ ' ' ^ lli=0 2r(2(m+j)) 

oT./. .N n"^'^ r(i+2j)^r(3+2i) 

_ 2r(4m - 4) llj=o 2r(2(m+i)-2) 



r(2m - l)2r(2m + 1) n™"^ r(i+2j) 
^ ^ ^ ^llj=o 2r(2( 



'r(3+2i) 



(2(m+j)) 

2r(4m - 4) "^^ r(2(m + j)) 

r(2m-l)2r(2m + l) -IJ- r(2(m + j) -2) 

2r(4m - 4) r(4m - 2) 

r(2m - l)2r(2m + 1) T{2m - 2) ' 
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To bound r(2m.-iV'r(2?n')r(2m.-2') f^°™ above we use the binomial relation (see e.g. [3, 
Eq. 6.1.21]) for the Gamma function 



^wj T{w + l)T{z-w + l) 
and the inequality (see e.g. [16, Eq. C.5]) 

{-r<(')<c-r (9.8) 

to obtain 

r(4m - 4)r(4m - 2) _ r(4m - 4) r(4m - 2) 

r(2m - l)2r(2m)r(2m - 2) 



r(2m - l)r(2m - 2) 


r(2m- 


- l)r(2m) 


/4m - 5\ /4m - 3\ 






V2m - 2j \2m - 2) 






/4m - 5\ /4m - 3\ 






V2m - 2) \2m - l) 






, Am - 5s2m-2/ 4m ■ 
^^2m-2'' ^"2m- 


— 3s 2m- 

-V 


-1 




(2e)2™-2(2e)'™-' 






(2e)4m-3 






\{2er-\ 







Recalling (9.6) establishes the assertion for n = 2m — 1. 

Now consider the even case n — 2m. We utilize the bijective principal representation 
(jy^i : A™ X T™' — ?■ Int{M'^'^) defined in (5.7) and, proceeding as in the odd case, we 
obtain 



^ ^ > l-2mr™"^^"^'^'^^ 



VolfM"^"^^ 5m(3, 1,2) 



Using Selberg's formulas (6.1) we compute 

5'm-l(3, 1,2) _ 11^=0 2r(2(m+j)) 



m-2 r(l+2j)r(3+2j)2 



5^(3,1,2) T-rm-l r(l+2j)r(3+2j)2 

^ ' ' ^ lli=0 2r(2(m+i)+2) 

■p-fm 

=0 2r(2(m+j)) 



oT./. o^ n"~^ r(i+2j)r(3+2j)^ 
2r(4m - 2) llj= 



r(2m - l)2r(2m + 1) n™"^ r(i+2j)r(3+2j)2 

^ ^ ^ ^■llj=0 2r(2(m+i)+2) 

2r(4m - 2) 7t r(2(m + j) + 2) 



1) -^i„ r( 



r(2m-l)2r(2m + l) JJ- r(2(m + j)) 

2r(4m - 2) r(4m) 

r(2m - l)2r(2m + 1) r(2m) ' 
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We will now use the Beta function 



B{a,b):= I t"-^ {I - tf-^ dt , a>0,6>0. 
Jo 



r(2m - l)2r(2m)2 B{2m - 1, 2m - 1) B{2m, 2m) 



< 



< 



2m — 1.^ n2m 



,4m,— 2 r)4m 



2m .^_o2m 



Am— 2 ^"''r)4m 



Finally, we apply the inequality 



< -m228™ 
- 16 



^2 < g/e 4„ 



'2' 
from Proposition 10.1 to conclude that 



2r(4m - 2)r(4m) 1 , g^ 

r(2m-l)2r(2m)2 - 8 

•C /'_\4mr)8rr 

= (26)'^" 



thus establishing the assertion for n = 2m. 



9.5 Proof of Lemma 4.4 



(9.9) 



and the identity 

see e.g. [', Pg. 258], where F is the gamma function. 

To bound r(2m-i"')^r(2?n)^ from above we use the inequality 

, , F(a)2 4 p„ 

from Proposition 10.2 to obtain 

F(4m - 2)F(4m) 1 1 



According to Chang, Kemperman, Studden [I 1, Thm. 1.3] one can show, using Skibin- 
sky's canonical coordinates for the moment problem [ ], that the uniform distribution 
on M" marginalizes to a Beta distribution corresponding to B{n, n) (see (9.9) and (9.10)) 
on the first moment. Consequently, 

yo/(gEM":gie [1-^,1]) _ 1 [' ,„_i., ,.„_i ,, 

Vol{Mn) - B{n,n)Ls ^' '> 

= h{n,n) 
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where Is{n,n) is the Incomplete Beta function (see e.g. [■>, Pg. 258]). Using the binomial 
relations (9.7) and and (9.8), for the upper bound we obtain 

nB{n,n) 
r(2n) 



= 5" 
< 5"ie 



nr(n)2 
2n-l 

n 
2n - 1 \ " 



n 
< 6"(2eV 



and for the lower bound 



''("■"' ■= i(;b) /'""'' -*>""* 



> TTT— ^(l-<5)"~' / t"~^dt 



B{n,n)^ Jo 

nB[n,n) 
2n —1 \ ,_ ^^„_l rn 



■;)(i-*)' 



n 

> (^^^)"~V-^)"-5" 

V n — 1 / 

> 2"~i(l-(5)"-^(5" 

> 5" 

where the assumption 5 < ^ was used in the last step. 

9.6 Proof of Proposition 5.5 

The following identity of Karlin and Shapley [ , Proof of Thm. 6.2] will be useful in all 
the Jacobian determinant calculations of this paper: For ti < si < ■ ■ ■ < tm < Sm, we 
have 

Qm 

— TT— A(tl, Si, ..,tm, Sm)\isi,..,s^)=iti,..,t^) = A^(t) . (9.11) 

Ubi ■ ■ ■ Ubm 

We can develop the upper and lower configurations simultaneously, by introducing a 
point to £ {0, 1} and representations (pto where when to = we have (J)q = 0^^ defined in 
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(5.15) and when to = 1 we have </>i = (/>q„ defined in (5.16). So, let us use this notation 
and a change of indices, and consider the two maps 

^™-i^/nt(M2™-i), to = 0,1 



^to 



A™ XT'' 



defined by 



m—l 



m—1 



<Pt,{\t) = ^(j^A.Jt^. +(l-^A,)<5i.) 

3=0 j=0 

m—l m—l 

j=0 j=0 

In this notation. Proposition 5.5 becomes 



2m- 1 



Proposition 9.1. For fo = 0, 1 we have 



m—l 



\det{d^t,)\{X,t) = Jt,{t) J] A, 



where 



m—l 



m—l 



Jtoit) =\to-t*\ii i^j - **)' n (*j- - *o)' ■ ^t-iit) 

3=1 3=1 

The differential of (j)to is determined by 



(9.12) 



''to 



dXi 



t]-tl, j = 0,..,m-l 



and 



dt 



*° -iXjff' i = l,..,m-l 



for i = 1, ..,2m — 1, from which we conclude that 

m—l 

\det{d(t>u,)\ = \Jto\ n^i 



where 



Jt 



to 



to — t^, 
'-0 f-* 



tl — t* 
t?-t? 



1 

2ti 



2m-l ,2m-l ,2m-l ,2m-l 



'0 



t: 



1 

«0 



1 

1 

h 



tr^^ {2m-l)ti 



1 

2ti 



2m~2 



2m-l +2m-l ,2m-l 



t, 







tf"-^ (2m - l)tf 



2m-2 



1 

^m— 1 
2 



t 



m—l 



m—l 



m—l 



,2m-l _ .2m-l 
''m-1 ''* 



1 

^'■m- 1 



(2m - l)C_f 




1 

•^'•m,— 1 



2m-2 



C-1^ (2m-l)C-l 
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To evaluate Jt^ for to = 0, 1, let si, .., Sm-i satisfy tj < Sj < tj+i, j = 1, .., m — 1 and 
define the Vandernionde determinant 



J{si,..,s^ 



m—lf 



1 



1 

to 



1 



Si 



j.2m-l j.2m-l .2m-l „2m-l 
'-* 1^0 1^1 b^ 



1 


„0 


tm-1 


Sm— 1 


f2 


q2 


.2m-l 


„2m-l 


''m-1 


*m-l 



and observe that the multilinearity of the determinant shows that 



Jt 



d 



\m—l 



to 



J\Si, .., Sm-l)|(si,..,s„_i)=(ti,..,t„_i) 



9si • • • dSm-1 
To evaluate this differentiation, observe that 

J{si,--^Sm-l) = ^{t*,to,ti,Si,...,tm-l,Sm-l) 

m— 1 m—1 

= (*o - **) ( n ^^1 - **)(^j - **)) ( n (*j - *o)(^j - *o) 

A(ti,Si, ...,tm-l,Sm-l) , 

from which we conclude that 

am—l 



Jt 



to 



9si • • • dSm-l 
m—1 



j{Sl, .., ■Sm-lj|(si,..,s„_i)=(ii,..,t„_i) 
m—1 



{to - i*) • n (*j - **)^ n (*j - *o)' 



9^ 



,m—l 



dsi--- dSm-l 

A(ti,Si, •.■,tm-l,'Sm-l)|(si,..,s™_i)=(ti,..,t„_i) 

Using the identity (9.11) we conclude that 

m—1 m—1 

jio = (io - i*) n (^j- - **)' n (*j - *o)'A^_i(t) 

thereby proving Proposition 9.1 and therefore Proposition 5.5. 



9.7 Proof of Proposition 5.6 

To simplify notation, let (j)i := (/>g; defined in (5.18) and (/)„ := (/)g„ defined in (5.19). We 
begin with the lower representation (pi. The differential of 0/ is determined by 

— -^ =f-—f i = \ m 
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and 



dti 



i>^jt] ^ j = l,..,m 



for i = 1, .., 2m, from which we conclude that 



where 



Jl 



\det{d4>i)\ = \Ji\llXj 
i=i 



tl 


-u 


1 


^m 


-u 




1 


tl 


-tl 


2ti 


•• tl 


-tl 


2tm 


+2m 
^1 


^2m 


2771*2™- 1 . 


^2m 


j.2m 


2mtl^-^ 


1 


1 





1 









t* 


h 


1 


*m 


1 






tl 


^1 


2ti 


• t2 


Ztffi 






t*™ 


j.2m 
^1 


2mt2™-i .. 


j-2m 


2mtl^ 


-1 





To evaluate J'l, let si, ..,Sm satisfy tj < Sj < tj+i,j = 1, ..,771 and define 



>-/ ("51, .., S-uiJ .- 



1 1 1 

t* tl Si 

fj. El Si 



2m j.2m ^2m 



1 


1 


^m 


^m 


t^ 
^m 


4 


1.2m 


2m 
^m 


^m 



^1 ^1 
and observe that the multilinearity of the determinant shows that 



Jl 



J [Sl, .., Sm)\(si,..,Sm) = {tl,..,tm) 



dsi ■ ■ ■ dSm 
To evaluate this differentiation, observe that 

«y (,Si, .., Sfn) — Ls[t^:, tl. Si, .., tfYi, Sm) 

and, using the recursion relation of the Vandermonde determinant, we obtain 

>-^ V"^!) ••) '^m j — ^it*, tl. Si, .., tm, S,nj 
m 

— II (J-j i^*)\Sj I* j ■ Ls[tl, Si, .., Em, Sfn) 
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from which we conclude that 



dsi ■ ■ ■ ds„ 

m 



o/ (Sl, .., Smj|(si,..,s„)=(ti, ..,*„) 



Q-n 



dsi ■ ■ ■ dsf 



-A{tl,Sl, ..,tm, Sm)\{si,..,Sm) = {ti,..,tm) 



i=i 
and therefore 

m 
3 = 1 

thus estabhshing the lower identity. 

Now, for the upper representation 0„ := (/)g„, the differential of (pu is determined by 



5An 



- = 4-tl, i = l,..,m-l 



d\, ^ 



dXr 



l-t 



and 



^ = ^W' ^• = i>-"^-i 



for i = 1, .., 2m, from which we conclude that 



m— 1 



\det{d(f>u)\ = [Jul n ^i 



where 



Uu 



-/■2 +2 _ J.2 



1 

2ti 



i^m—l 



t 



m—1 



_^2m i2m_^2m 2mt?"-^ 



1 
•^I-m— 1 



l-t* 
l-t? 



j.2m j.2m n j.2m— 1 i j.2m 



t* tl 



1 
2ti 



m—1 



^2m ^2m 2mt^ 



2m- 1 



m,— 1 



2t 



m—1 



2m- 1 



'-m-l ^'"-''m-l ^ 
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To evaluate Ju-, let si, .., s^-i satisfy tj < Sj < t,+i, j = 1, .., m — 1 and define 



J{si,..,Sm-l) ■■-- 



t* tl Si 

f2 f2 2 

''* '■1 *1 



i2m j.2m „2m 
f-* ''1 *! 



'-m— 1 Sj^j—i i 



m—1 ^m—l 



j.2m „2m i 

''m-1 ''m-l -"^ 



and observe that the multilinearity of the determinant shows that 

am— 1 



•Ju 



dsi--- dSm-l 
To evaluate this differentiation, observe that 



J{Sl, .., •Sm-l)|(si,..,s„_i)=(ii,..,tm_i) 



J{si,..,s 



m—1 1 



m—1 m—1 

t* n *^ n 





1 




1 




1 




1 




1 




1 


1 


1 




t* 




h 




Si 




t2 




S2 




tm-1 


Sm— 1 


1 


SJ 






























j.2m- 


-1 


.2m- 
^1 


-1 


2m- 
^1 


-1 


,2m- 


1 


,2m- 
•'2 


1 


.2m-l 
■ "-m-l 


2m-l 
*»n-l 


1 



That is, we have 



m—1 m—1 



J{si, ..,Sm-l) — t* _[_[ tj _[_[ Sj ■ A(t*,ti,Si, ..,tm-l,Sm-l, 1) • 

We use the recursion relations 

m—1 
A(t*,tl,Si,..,tm-l,Sm_l,l) = (1 - t*) Y\^ (1 - tj)(l - Sj) • A{t^,ti, Si, ..,tm-l, Sm-l] 



j=l 



and 



m—1 



A(i*,il,Si,..,tm-l,Sm_i) = JJ (tj - t*)(sj - t*) • A(ti,Si,..,t„_i,Sm-l) 

to obtain 

m—1 m—1 

J{si,..,Sm-l) = t*(l-t*) JJ tj(l - ij)Sj(l - Sj)- JJ (tj - t*)(Sj - t*)-A(ti,Si,..,tm-l,Sm-l) 

Consequently, the identity (9.11) implies 

am- 1 



Ou 



J [Si, ■■, Sm-lj|(si,..,s™_i)=(ti,..,tm_i) 

m—1 ry^n—l 



dsi--- dSjn-l 

m—1 

i*(i - i*) n *i(i - *i)^ n ^*j ~ * 



d" 



dsi--- dsm-i 

A(t*,ti,Si,..,tm-l,Sm-l)|{si,..,s™_i)=(ii,..,t„_i) 



m—1 



m—1 



t*(i - u) n (h - 1*? n *'(! - h)' ■ ^'m-iit) 

j=i i=i 
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establishing the upper identity and thus completing the proof. 

9.8 Proof of Proposition 6.1 

For the first assertion, let n = 2m — 1. As in the proof of Proposition 5.5 we find it 
convenient to analyze the upper and lower configurations simultaneously, by introducing 
a point to G {0,1} and the volume filling representations (j)to,to G {0,1} where when 
to = we have (po = (p'^i defined in (5.15) and when to = 1 we have ^i = cj)^^ defined in 
(5.16). In this notation, from (9.12) we have 

m—l m—1 

<Ptoi^O,--,^m-l]h,--,tm-l) = / , ^j't) + (^ ~ / ^ -^i)^* 

i=0 i=0 

m—l 

j=0 

m—l 

= Ao (tj, -t:) + ^ A, (tj-t:)+t: (9.13) 

for to = 0, 1 and Proposition 9.1 expresses the Jacobian determinants as 

m—l 

\detidcf>t,)\{X,t) = Jt,it) J] A, 

i=i 

where 

m—l m,— 1 

jtoit) = iio - t*i n (*j- - **)' n (*^- - *o)' ■ ^--iw • (9-14) 

i=i i=i 

In this notation, the modified change of variables formula (5.28) becomes 



Af2m-1 tn = l 2'^^'"^^'""^ 



for i > 0. Therefore, we conclude that 



/ 



M2™-i to=l 2'^"^'"^^™'"^ 



1=1 2-^^'"xT'™'~^ 

„ m—l 

3 = 1,2"^^^ ^^ j = l 

„ „ m—l 

^2-''^ •^^™ f=l 
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(9.15) 



(9.16) 



Performing the A*" integration, using the identities /^^ >\n-i HiLi ^i^^A = (2^7 /a™ ni^i ^i^^^ 
(2^, and /^,„ UT=i^ ^id^ = (2;;^' we obtain 

„ m— 1 „ m— 1 m— 1 

/ < n ^i = / (^0 {th - ii) + E ^i (*^- - **) + **) n ^^-^^ 



m— 1 



1 • 2 ' - 1 . 

(2;^ ^*° " **) + (2;;^ ^ ^*^' " **) + (2m -1)!** 



1 . 2 "-' 



rtn + 



V f • + — 



(2m)! '^ (2m)! ^ ^ (2m)! 
Consequently, for i > 1, we have 



^0 = 1,2'^-' -^^ 7 = 1 



/\/f2Tn — 1 






m— 1 



(2m)! 
and, for i = 



L-, ^' + p^'- X»-. (■^» + -^0 + p^ X,„-. (E *5) (■7» + •^' 



« = » ' iM / (Jo + Jl) (9.17) 

4^2m-l (2m — Ij! Jj-m-l 



which we aheady knew from (6.7). Combining the two, we obtain for i > 



JV/2m-l {2m)\ Jj'm-l {2m)\ Jrpm-1 (2m)! * Jj^m-l 

„ » m— 1 

+ (2;^i,...(g'5)(-7» + -7') ("«> 

and the substitution of the volume equality (9.17) (that is, (6.7)) yields the assertion in 
the odd case. 

For the even case, let n = 2m, and let us simplify notation by denoting the volume 
filling representations by cpi := (j)'^i and 02 := (peu defined in (5.18) and (5.19) so that, in 
this notation, 

is defined by 



m m „ 

2m 



(t)l{\u..,\m-M,..,tm) = (E^i4 + (1-E^i)**)iri 
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and 

by 

m—1 m „ 

(p2{^0, ■■,^m',ti,..,tm-l) = (/_^^jt] + ^m + {^ — /_^^j)tl],_ ■ 

j=l j=0 

From Proposition 5.6 we have 

m 

|det(#i)(A,t)| = Ji{t)llXj 

i=i 

771—1 

|(iet(#2)(A,t)| = J2{t)'[{Xj 

i=i 

where 

m 

^iW = Uitj-Uf-Aiit) 

m— 1 m—1 

In this notation, the modified change of variable formula (5.28) becomes 

/ qi= f 4\dM+ I 4l#2|. (9.19) 

We evaluate the two integrals in (9.19) by 

„ . m 

/ 'All^ll = / <l^l([\^,)Jl 

JA"^xT"^ JA'^xT"' -^ 

« „ m 

= / (/ 0invA)^i 

and 

» » m—1 

/ 4\dcpi\ = / 'A2(n^^)^2 

JA'"+ixr™-i Ja^+IxT™-! 1 

„ „ m—1 

= / (/ '^invA)^2 



i=i 
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Performing the A'"^^ and A*^ integrations, using the identities J^™. Ai H^i '^jdX = 

(2m+iy. ' Ja™ '^m-2 llj=l ^i"'^ = (2m-l)! ' Ja™ 11j = 1 '^i"^ = (2m-2)! ' Ja™ "^m-l lli=l ^i"-^ 

(2^' /a'" n^i ^id>^ = (2^' and /^„ llT=i^ ^id^ = (2m-i)! ' we obtain for i > 1 

„ m „ m mm 

/ c/^iHAjdA = / (^A,4 + (l-^A,K)nA,rfA 

„ m m 

(2m + l)!^^^ *' {2m)\* 



J 

m 



1)!^ ^^ (2m + 1) 



(2m + 1)! ^ ^ (2m + 1) 



-tl. 



and 

m—l „ m— 1 m m—1 



n I I li J- p IIV J. I I It IIV J. 

/ </.^ n A,dA = / (V A,4 + A„, + (1 - V \,)t) n ^i^^ 

■/a™+i J=i -^^'"^^ Stt jt^ ^ J=l 

„ m—l m—l 

= / ( E ^i(*J - **) + ^™(1 - 4) + (1 - Ao)t:) n ^i^^ 

2 m.—l ^ ^ ^ 

(2m + l)!^^^ *^ (2m + l)!^ *^^(2m)! (2m + l)!- 



m—l 



_ V t^- + 3^ + i f 

1)! ^ ^ (2m + 1)! (2m + 1)! *' 



(2m + l)!^^ (2m + 1)! (2m + 1)! 

For i = 0, (9.19) imphes 

Vol{M^^) = -^ / Ji + -^ / J2 
(2m)! j-rm (2m)! Jx^n-i 
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so that for i > 1 we have 



4>\\d4>i\ + I (P2\(^(P2\ 



ir- V (2m + 1) ! ^ *^' ^ (2m + 1) ! * V ^^ 
^ ir— 1 V(2m + 1)! ^ *^' ^ (2m + 1)! ^ (2m + l)!*V ^^ 
= ir™((2m + l)!^*^)^^ 

. „ m-l ^ 

^ ir^-i l(2m + l)! ^ *^' ^ (2m + 1)! J ^^ 



+ ** VoliM^"") 

2m + 1 ^ ^ 

so that for z > 1 we conclude 

r2m 



m— 1 



JA/2m 2m + 1 

Combining with the result (9.20) for i = we conclude 

JM^m 2m + 1 

f^ p TTl or 171— L ^ (^ / .\ ^ 

(2m + 1)! ir,,. ^ ^^'^^ ^ (2m + 1)! /r— i ^ ^^'"^^ ^ (2m + 1)! 7^—1 '^^ 
establishing the assertion in the even case. 

9.9 Proof of Theorem 7.1 

Recall the identity J^oi{0,t) = 0. Then subtracting the volume identity (6.7) 

from itself evaluated at i* = 0, we conclude that 

/ {Joi(.t.,t) + JUt*,t) - Jou(.0,t))dt =0, (9.21) 
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that is, 



/ n{U,t)dt ={). 

J lm~\ 



We now do the same subtraction for all the moments. To that end, recall the convention 
0*^ = 1, and observe that, for i > 0, the identity 

= (2^).;1?- 1)1 i- '^°'('""* ^2.,)l(m-l)! /„-. ■^»('-"* 
from Proposition 6.1, evaluated at t^, = becomes 

T2m-1\ 



[ ,. - '-^VoliM^ 



^'^'^ ^ ^o1(o, t)dt + ,,^,„i ,„ / J,^(0, t)dt 



+ 



(2m)!(m-l)!yj™-i °'' ' ' {2m)\{m - l)\ J jr^-i 

2 



-^ y^^_^ st^ ( j,^; (0, t) + j,^„(o, t)) dt . 



(2rn.)!(m 
Subtracting from (9.22), using the identity J7o'](0,t) = 0, we obtain 

2m 2m 

and applying the volume identity (9.20) we obtain 

-^yo/(M2™-i) 

2m ^ ^ 



Jim.-i (2m)!(m - 1)! Jjm-i 



(2m)!(m- 1)! 

— — -^ — / T.fn{U,t)dt 

(2m)!(m — 1)! J/m-i 



and the subtracted volume identity (9.21) we obtain with a change of sign 
'/ol{M'^'"''-^) 



-^Vol{M'^"'-^) (9.23) 



(2m)!(m- 1)! 



/ J^{U,t)dt + ——^ — / J,'i{t,,t)dt 

Jjm.-i (2m)!(m - 1)! Jjm-i 



(2m)!(m- 1)! jjra-i 



+ TKZXTTZ: TT7 / ^fnit^,t)dt. 
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Then, if we let (j){s) := X^j^ </'iS* be a polynomial of degree n = 2m — 1, summing 



over each identity in (9.23), we conclude that 
Wol{M 



'^(**)T7„U,^2m-l^ 



2m 

2m- 1 



TT7 / {Joi it* , t) + J^^ (t. ,t))dt+Yl ^i TT^ TTT / ^oi (t* , t)dt 

m. — 1)\ J jm~i ' ^ [2m,y.[m — 1)\ Jjm-1 



- '"^2m)!( 

+ TT^ nT / {m{t)nt*,t))dt 



1 f ^""""^ 1 f 

[2m)l[m, — 1)\ J jm-i ^ (2m)!(m — Ij! jjm-i 

+ T^T^TT^ nT / {m{t)nt*,t)dt. 

[2m.)\(m — 1)\ Jjm^i 

Since 0(0) = (/>o and (/>(!) = X^jl^" i;^j the assertion follows by multiplication by 2m,. 
The even case proceeds in the same way, but since it is a little different we have included 
it in Section 10.2 in the Appendix. 

9.10 Proof of Theorem 8.1 

For the first assertion, let n = 2m, — 1 and consider the integral formula (9.20) 

V'-m'"'-') = (2„.-l)!(„.-l)! /._, (•7»'('- ') + -^^C- '))"* 
for the volume in terms of the canonical representations. From the definitions 

ra—l m— 1 

rra— 1 m— 1 

JLit*,t) = {l-u)ll{t,-ufll{l-t,)^.Ai_,{t) 
of Proposition 5.5, we obtain 

j=i 

^J„'L(t.,t)k=o = -{l + 2^t-')llt]{l-t,fAt^_,{t). 

j=i 

Differentiating the volume formula with respect to t* at t* = 0, we obtain 

» = (2„. - i)K,. - 1). /.,- (l:-^.sfe.')i'-° + |:^^('..')i..-o)^t 
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and therefore 

m— 1 „ m—1 



p lib J. n Illy _L 

I Tm—1 I Tm — l 

Tl— 1 

/m— 1 
,=1 



m-lit)dt 



m—1 



Tni — l 

^ i=i 



/m— 1 
St-1 • n t]{l - tjfAt-i{t)dt = 5^_i(5, 1, 2) - S„_i(3, 3, 2) 

.7 = 1 



m— 1 

+ 2 ' 
from which we conclude that 

m— 1 

2 

Changing m, i-^ m + 1 finishes the proof of the first assertion. 

For the second assertion, let n = 2m and consider the integral formula (6.8) 

VoliM^"") = . } , / J::i{u,t)dt+ . .,/ — / J,liu,t)dt. 

{2m)\m\Jim {2m)l{m - 1)\ J jm-i 

From the definitions 

m 

J=l 

m—1 m—1 

i=i j=i 

of Proposition 5.6, we obtain 

8 "^ 

^Jel(**,*)lt,=o = -2J:t-'l[tlAUt) 

i=i 

rs m—1 

i=i 
Differentiating the volume formula with respect to t^^ at t* = 0, we obtain 

= ^9 M I / 7^^e1(**,*)lt.=0rft + 7:^^77 7T7 / ^ Jen(**, *)lt.=odt 

(2m)!m! jjm at* (2m)!(m — Ij! J/m-i ai* 

and therefore we conclude 

/m „ m—1 

j:t-'-lltlAtit)dt = m l[t^{l-t,)^-At,_,it)dt 

= mSm-i{5,3,2), 
finishing the proof of the second assertion. 
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9.11 Proof of Theorem 8.2 

First note that the definition (7.1) 

and the definitions of J'^i and J'^^ from Proposition 5.5 imply that 'H{-,t) G n^™~^, t G 
/™-i. Therefore, it follows from (7.2) that 7i(-,t) e n^^-^, t G I™-^ Now, it follows 
from Theorem 7.1 that 

0(t*) = / {^(t>){t)n{u,t)dt, ^eul""-^ 

Jjm-l 

which expanded becomes 

2m-2 „ 

4>{u) = Vpi(t*)/ {Y.4>){t)hj{t)dt, 0Gn2"-S 

in particular, by choosing (j) := pj., k = 1, .., 2m — 2, 

2m-2 .. 

Pk{t*) = V" Pi(i*) / {T.pk){t)hj{t)dt, k = l,..,2m-2, 

from which we conclude 

T,pj ■ hk = 5jk, j, k = l, .., 2m - 2 
-1 

establishing the assertion. Furthermore, from this and the symmetry of /ifc, , /c = 1, .., 2m— 
2 with respect to the action of the symmetric group, we also conclude 

(m - 1) / pj ■ hk = Sjk, j, k = l, .., 2m - 2 

establishing Corollary 8.3. 

9.12 Proof of Theorem 8.4 

From the orthogonality relation (see e.g. [5, Eq. 12.110]) 

J,^^^'' = 2k + l{k-2)?^' 
and the definition 

■ yo?(M2'»-i) (2m - l)!(m - 1)! 
of the scaling of the kernel (7.1) 

m*,t) = JL{t) - J^AU.t) - Jolit.,t) , (9.24) 
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we can compute the coefficients /i^ in tlie expansion 

2m-l 

n{U,t) = Y^ hj{t)Qj{Q, {U,t) e / X /"-I 

as 

, . X 1 2 (2k + l)(k-2)\ r , ,^ , , 

(9.25) 
and then apply Theorem 8.2 to obtain the assertion. To that end, for k = 2, .., 2m — 1, 
to compute 

ln;t)Qk{-) 

we use the decomposition of (9.24) of 7^ and compute the values Jj JSu,{t)Qk[')-, fj 'Joi('^^)Qki'), 
and fj J^ou{'^'^)Qk{-) separately. For the first term, observe that /j Qfe = 1 + (—1)^, k = 
2, ..., 2m — 1 from Proposition 9.2, so that 

j^JL{t)Qk{r)dr = JL{t) JQk{r)dr = (1 + (-1)^) J^Jt) . (9.26) 

For the second, fjJ^oii-,t)Qk{-), we expand J'^i{-,t), defined in (5.5), as a polynomial 
for fixed t and then utilize the values of the integrals Jjr^Qki'>')dr for the monomials 
r^ ,j = 1, .., 2m. — 1. To that end, define 

m—l 

J:ii;t):=t*l[{t,-Uf (9.27) 

i=i 
so that 

m— 1 

Then from the definition of J"^ of (5.4) we have ^^(t) = UfJi^ t] ' ^m-i(*) so that 

Joli;t) = J^l{;t)X'l{t). (9.28) 

The generating function identity for the elementary symmetric functions Cj is 

m— 1 m—l 

i=l j=0 

and squaring it we obtain 

m—l m—l 

llil + sUf = (^e,(t). 
i=l i=o 

2m-2 

= ^ e,(t,t)s^ (9.29) 

i=o 
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is-' 



2 



Therefore, by changing s i— > —t^ ^ we conclude that 

m—l 2m— 2 

Y[{tj-Uf = Y, ''^rn-2-j{t,t){-iyti (9.30) 

j=l j=0 

and therefore 

m— 1 2m— 2 

U n {t,-Uf = J2 e2m-2-j{t,t){-iyti+' 
j=l j=0 

and, relabehng, we conclude 

2m- 1 

Jo'iit*,t) = - Yl djm-iyti (9.31) 

i=o 
where 

dj{t) := e2m-i-j{t,t) = Yl ej,{t)ej,{t), j = I,. .,2m - 1 (9.32) 

jl+J2=2m-l-j 

and do(^) '■= 0. Note that dj is a symmetric polynomial of degree 2m — 1 — j. 

In particular, from the definition (9.27) and its resulting polynomial expansion (9.31) 
we have 

2m— 1 m—l 

J5(l,t) = - Y Mt)i-^y = n i^-h)'- (9-33) 

j=0 j=l 

The following proposition computes the values of the integrals of the Legendre polyno- 
mials against the monomials, and we observe that Jjr^Qkif)dr = 1 for 1 < j < k and 
// r^Qk{f)dr = 1 plus a term when j > k. 

Proposition 9.2. For k = 2, .., 2rn- — 1 we have 

j Qf,{r)dr = l + {-lf 

fr^Qkir)dr = l, l<j<k 

f 3n ( w 1 (j + fc + F)r(j + 2)r(j) . 

/ r-'Qkirjdr = 1 — — ; — — ; -, i > k . 

Therefore, since do = 0, we obtain from the polynomial expansion (9.31), Proposition 
9.2, and (9.33) that 

/2m— 1 2m— 1 

J^i{r,t)Qk{r)dr = ^ (-l)^a,,d,(t) - ^ (-l)^d,(t) 
j=k j=0 

2m- 1 

= Y(-^y^Jkdj{t) + X%ht) (9.34) 

j=k 
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where we recall the definition 

{j + k + A;2)r(i + 2)r(j) 



O-jk 



■ r(j + k + 2)r{j - A: + 1) ' 



and note that the lower limit in the summation is k and consequently, the first term 
is a symmetric polynomial of degree 2m — 1 — k. Therefore, multiplying by J^iit) = 
WT=i t^j ■ ^m-i(^) ^'^'i using the definition (9.28) and the identity 

m—l m— 1 

= JLit) (9.35) 

from the definition J^{t) = Hfci ^?(1 ~ ^jY ' ^m-i(*) from Proposition (5.3), we con- 
clude that 

„ 2m-l 

J^i{r,t)Qu{r)dr = J-P(t) ^ (-l)%^.,d,(t) + J-<fJt) . 
Let us designate the negative of the first term 

2m- 1 

kit) := J!i{t) Y, {-iy^'ajkd,{t) (9.36) 

j=k 

SO that 

j Joi{r,t)Qk{r)dr = -hk{t) + J^^{t). (9.37) 

For the third term J^Jou{'^t)Qj{-), we utilize the reflection symmetry 

of (5.17) and the reflection symmetry (see e.g. [ , Eq. 12.97]) 

Q,,{l-r) = {-lfQk{r) 

of the associate Legendre polynomials to compute the integral fjJ^oui'^t)Qji') i^ terms 
of//^ol(-.OQi(-)- That is, 

j JL{r,t)Qk{r)dr = j J^{1 - r,t)Qu{l - r)dr 

= Jj,Ur,l-t)Qk{l-r)dr 

= (-1)' [jo'i{r,l-t)Q,{r)dr 
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and therefore 

J Jouir,t)Qk{r)dr = (-1)^ J J^^(r,l - t)Qk{r)dr 

so that we conclude 

Jouir, t)Qk{r)dr = -{-l)'k{l - t) + (-1)' J'.TJl - t) . (9.38) 



ou\ 

I 



Putting all three terms together using (9.24) and the identities (9.26), (9.37), and (9.38), 
along with the symmetry 

we conclude that 



'Hir,t)Qkir)dr = / jP^{t)Qkir)dr - / J-{r,t)Qkir)dr - / Jol{r,t)Qkir)dr 
I Jl Jl Jl 

= (1 + (-1)'=) J?Jt) + h,{t) - Ji^(t) + {-\fk{\ - 1) - {-\fJL{^ - 1) 

= hkit) + {-!)'' ka-t). 

To finish, consider the functions hk{t) := hj^{t) + {—l)^hk{l — t). It follows from 
(9.25) that the basis coefficients h^ satisfy 

1 2 {2k + l){k - 2)1 ^ 

'' ~ Vol{M^"^-^) (2m - l)!(m - 1)! (fc + 2)! '=' ^-2,",2m-l. 

Moreover, Theorem 8.2 implies that {SQj, j = 2, .., 2m — 1} and {/i^, k = 2, .., 2m. — 1} 
are an L^(/™'~^) biorthogonal system. It therefore follows that 

^^ ^ hkT.Q, = yo/(M2™-i)(2m - l)!(m - l)!^^^i^-t^L_<5,,, j. A: = 2, .., 2m - 1 . 

Moreover, from the symmetry Qk{^ — r) = ( — 1) Qfc(r) and a change of variables we 
obtain 

/ hk{t)T.Qj{t)dt = / hk{t)T.Qj{t)dt + {-lf hk{l - t)^Qj{t)dt 

= I hk{t)T.Qj{t)dt + i-l)^ f hk{t)^Qj{l-t)dt 
= I hk{t)^Qj{t)dt + {-iy+'' [ hk{t)^Qj{t)dt 

Jjm-l Jlm-1 

= (1 + (-1)^+'=) I hk{t)T.Q,{t)dt. 

Since {j + k) mod 2 = [j — k) mod 2, the assertion then follows from the definition (9.36) 
of hk, the identity J^i{t) = Y\T=i t] • A^_;^(t), and the definition (9.32) of d^- Moreover, 
we see that for (j — k) mod 2 = 1 the vanishing of this integral does not depend on the 
function h^ but is instead a consequence only of the relative parity between Qj and h^ 
with respect to the operation of reflection. 
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9.13 Proof of Proposition 9.2 

We abuse notation by letting P^ and Qk denote the Legendre polynomials on the stan- 
dard set [—1, 1]. At the end we will change back to the interval /. We use [ :"', Eq. 7.127, 
pg. 771] 

(1 + xyPk{x)dx = — Z, , 7, J > 

1 > ^"^ > r(j + A; + 2)r(j - A; + 1) ' '' - 

and the definition Qfc(x) := (1 — x^)P^'(x) and integration by parts. Because of the poles 
of the Gamma function at and the negative integers, we conclude that 



/ {l + xyPk{x)dx = Q, j<k. 



Consider the function </>(x) := (1 + x)-' (1 — x^) = 2(1 + xy^^ — (1 + xy~^'^, which has 
the derivatives cp'ix) = 2{j + 1)(1 + xy - [j + 2)(1 + xy+^ and ^"(x) = 2{j + l)j{l + 



v-i 



(j + 2){j + 1)(1 + xy. Since, for j > 1, we have <I)'{-1) = and (p'il) 



-2J 



+1 



and Pfc(l) = 1) we obtain 



[l + xyQkix)dx 



{i + xy{i-x')pl:{x)dx 



(j){x)P'^{x)dx 



<PPL 



(j)' {x)Pl{x)dx 



-1 



(t)'{x)Pl.{x)dx 



-1 



b'Pk 



-1 



+ / (f)" {x)Pk{x)dx 



2^'+^+ / (j)" {x)Pk{x)dx 

2^+1 + 2(j + l)i / (1 + xy-^Pk{x)da 

-{j + 2){j + l)j {l + xyPk{x)dx 



from which we conclude that 



/ {I + xy Qk{x)dx = 2^+^ , l<j<k. 
For the case j = 0, defining </>(x) := 1 — x^, we instead have </>'(— 1) = 2, (j)'{l) = —2, 
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Pfe(l) = 1, and Pfc(-l) = (-!)''• Using /^^ Pk = 0,k>l, we obt 



am 



1 ^1 

Qk{x)dx = / {I - x^)P'^{x)dx 



-1 



-1 

1 



(j){x)Pl!{x)dx 
(j)' {x)Pl{x)dx 



1 fi 

-0'Pfc + / (/)"(x)Pfc(x)dx 

-1 J-1 



-4>'Pk 



-1 



-2Pfc(l)-2Pfc(-r 
2(1 + (-1)'). 



On the other hand, for j > k we have 
-1 



{l + xyQk{x)dx-2^ 



1+1 



2(i + l)j / {l + xy-^Pk{x)dx-iJ + 2){j + l) f {l + xyPkix)dx 



2^+ii(j + i)r=^(i)' ^ 



2^ 



r(j + fc + i)r(i - A:) 
jij + i)r2(j) 



2^+^(i + 2)(j + l] 
j(i + 2); 



r(i + fc + 2)r(i - fc + 1) 



Tij + k + 2)r(j - A; + 1 

^,+i i(j + i)(i + fc + fe')r'(i) 
T{j + k + 2)r(j - fe + 1) 

,,+i (j + fc + fc^)r(j + 2)r(j) 

T{j + k + 2)T{j -k + 1) 



T{j + k + 2)r(j - A: + 1) 
iJ + k + l)iJ-k)-JiJ + 2) 



and therefore 



Translating to the unit interval with the map / i— > [—1,1] defined by x = 2r — 1 we 
obtain the assertion. 
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10 Appendix 

10.1 Proof of even case of Lemma 4.3 

We utilize the bijective principal representation 0^^ : A™ x T™ -^ Int{M'^"^) defined in 
(5.7) and 

m 

\det{dreum,t) = jiMll^^ 

where 

m 

from Proposition 5.4. 
Fix t^, G (0, 1) and let 

It follows that 

Ml"" D(/)^;(A'" xTf) 

and therefore 

VolUll"') > yo/(0P (A™ X T^)) . (10.1) 



Using the identity J^„ TVILi ^i^X = /2m)\ ^ ^^ compute the righthand side using the 
change of variables formula as 






(2m)! y^-'' 



1 
{2m)\m\ J J 






where 

ir := {(ti, ..,U G I*" : t, i Bs,j = 1, ..,m}. 

To bound this from below we bound the integral over (/^) from above. To that end, 
let 

/,- := {(ti, ..,tm) e /" : tj E Bs}, j = 1, ..,m, 

so that 
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Therefore, using a union bound, we have 

„ ,. rn 

JilY') J (it) 7-1 



^s ) j=l 



and so obtain 



lltlAt^{t)dt 

I iTn Ttrt 



m 
III 



m 

lltlAtXt)dt 

„ m 

= rn T\t] T\ {tk- tj)Uh ■■■dtr, 

''^S,ij=l l<j<k<m 
,. rn 

< rn Wt] W {tk-tjYdti---dtr, 

'^h,ij=2 2<j<k<m 
„ m 

Vol{Bs) T\t] ■ Ai_^{t)dt 



m 

'^-' j=2 

= mVol{Bs)S„,^i{3,l,2) 
< 2mSSm^ii3,l,2) 



Vol ((^eilA"^ X T™)) //" UT=1 ^ ■ ^fnit)dt 

-5'm-l(3, 1,2) 



> 1 - 2m6- 



^rny'J, i, 2) 
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Using Selberg's formulas (6.1) we compute 

T-rm-2 r(i+: 

•Sm-llS, 1,2) _ llj=0 2r(2(m+j)) 



m-2 r(l+2j)r(3+2jf 



5^(3,1,2) Y\T=o 



2T{2{m+j)+2) 

-r\m-i r(i4 

-(2(m+j)) 



2T{4m-2) llj=o ~^r( 



r(2m - l)2r(2m + 1) n™"^ r(i+2j)r(3+2j)2 

^ ^ ^ ^J-lj=0 2r(2(m+i)+2) 

2r(4m-2) 7t r(2(m + j) + 2) 



l)2r(2m + l) J-J- 



r(2m-l)2r(2m+l) -fj- r(2(m+j)) 

2r(4m - 2) r(4m) 

r(2m - l)2r(2m + 1) T{2m) ' 

10.2 Even case of Proof of Theorem 7.1 

For the even case, n = 2m, recall the identity v7eu(0, i) = 0, and the volume identity 
(6.8) 

Vol^M^"') = —\^ f XiiU,t)dt + -—^ / Je'uit*,t)dt (10.2) 

[2m)\m\ Jjm [2m)\[m — Ij! J/m-i 



Now observe that, for i > 0, the identity 
Jm'2"^ 2m + 1 



^2m^ 



'°^^^ + ' ' ^e^Jt.,t)d. 



(2m + l)!(m-l)! 7/m-i 
2 



— — / sfj;^^(t*,t)(it + - — — ^- — / ^fj^^{u,t)dt 

ly.mljjm {2m + l)\{m - 1)\ J j,n-i 



{2m + 
from Proposition 6.1, evaluated at t* = becomes 

Jj^pm 2m + 1 

6o{i) + 1 



r2m\ 



(2m + l)!(m- 1)! y/m-i 



(2m + l)!m! J/m ""' (2m + l)!(m - 1)! J^m-i 

where we note that T,f = Y^T=i t!- in the integral over /*" and St' = Y^^Si t* in the 
integral over I"^~^. 
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Subtracting the two, using the identity J^ui^i't) = 0, we obtain 

Jeuit*,t)dt 



2m + 1 ^ ' 2m + 1 

6oii^ + 1 



(2m + l)!(m-l)! Jjm-i 
and applying the volume identity (10.2) we obtain 

ti 



2m +1 



i — VoliM^"") 






Then, if we let (j)(s) := Yli=o *^«** be a polynomial of degree n = 2m, summing over 
each identity in (10.3), we conclude that 



2m + 1 ^ ^ 

1 



,(2m + l)!(m-l)! Jj,^-i 

2m ^ , 



J.'„(«.,tMt- p,„^\),„, /_J/,(«..«)* 






Since ^(0) = (po and i;^(l) = X^j^g (pi the assertion follows by multiplication by 2m+l. 

10.3 Assorted Technical Results 
Proposition 10.1. We have 

m2<8(-)^", m>l. 
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Proof. We proceed by induction. The inequality is clearly true for m = 1,2. Therefore, 
suppose that it is true for some m >2. Then, since 

m + 1 3 /Cs 2 

< - < (-) , m> 2 

m - 2 - ^2^* ' 

follows from 6 < e^, we conclude that 



< 



{m+ if 



\Am 



< l^\'^Qf^\'im 



< 



m? 2' 

'£\4m+4 

thus establishing the inequality for m + 1 and finishing the proof. D 

Proposition 10.2. We have 

5(a,a) > -2~2", a>l 
a 

Proof. We have from the integral formula 

B{a,b) = I e-\l-tf-^dt 
Jo 

Jo 

> 2{^y~' Tf^-^dt 
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